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Abstract 


This  final  technical  report  summarizes  the  research  carried  out  by  our  group  concerned  with 
arcjet  plasma  propulsion.  Section  I  describes  the  non-equilibrium  plasma  model  and  a  detailed 
formulation  for  kinetic  non-equilibrium.  Section  n  describes  experimental  measurements  in  which 
measurements  of  electron  number  density,  electron  temperature,  ion  speed  ratio  and  gas  velocity 
are  carried  out  by  advanced  Langmuir  probes.  Section  m  outlines  our  research  with  a  modified  1- 
2  kW  class  arcjet,  where  surface-mounted  Langmuir  probes  are  used  to  investigate  the  constrictor 
and  nozzle  current  attachment  regions. 
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Overview 


A 


A  general,  two-dimensional,  steady,  non-LTE  (local  thermodynamic  equilibrium)  model  is 
presented  for  nitrogen/hydrogen  arcjet  thrusters.  Both  thermal  and  chemical  nonequilibrium 
processes  are  incorporated  to  describe  a  multi-temperature,  multi-species  N2/H2  plasma.  The 
current  and  anode  temperature  distributions  are  fully  coupled  to  the  plasma  flowfield  with 
realistic  boundary  conditions.  Thermal  nonequilibrium,  chemical  equilibrium  results  are 
presented  for  simulated  hydrazine.  Results  for  both  thermal  and  chemical  nonequilibrium 
hydrazine  and  hydrogen  arcjets  show  significant  thermal  nonequilibrium  near  the  anode,  and 
high  frozen  flow  losses.  It  is  concluded  that  nonequilibrium  chemistry  models  are  necessary  to 
capture  accurately  the  interior  distributions  of  the  plasma  flowfield. 

With  recent  advances  made  in  low  power  (1-2  kW)  arcjet  numerical  modeling,  validation 
of  numerical  models  through  experimentation  has  become  increasingly  important.  We  describe 
multiple  electrostatic  probe  techniques  for  characterizing  flow  conditions  at  the  exit  plane  of  a  1 
kW-class  hydrazine  thruster.  An  improved  electrostatic  quadruple  probe  technique,  which  is 
also  modified  to  account  for  plasma  gradients  over  the  finite  probe  volume  and  multiple  ion 
species,  provides  measurements  of  ne,  Te  and  uj/cm  at  the  thruster  exit.  Quadruple  probe  results 
are  also  presented  for  measurements  at  several  centerline  locations  immediately  downstream  of 
the  thruster  exit  plane.  A  spatially-resolved  time-of-fiight  electrostatic  probe  technique  is 
employed  for  measurements  of  radial  profiles  of  the  plasma  axial  velocity  ui.  The  quadruple 
probe  results  (Te,  ne,  ui/cm),  coupled  with  independent  measurements  of  plasma  velocity,  are 
used  to  estimate  ion  temperature  Ti=Tg  and  the  extent  of  thermal  non-equilibrium  at  the  thruster 
exit.  Exit  plane  experimental  data  for  ui,  Te  and  Tg  are  compared  with  computational  arcjet 
model  predictions,  showing  agreement  for  uj  and  Tg.  The  model  underpredicts  Te  and  a  direct 
comparison  of  ne  results  is  precluded  by  large  axial  gradients  at  the  thruster  exit. 

A  1 .4  kW  arcjet  thruster  has  been  designed  for  internal  probing  of  the  nozzle.  The  arcjet 
employs  a  single-piece  nozzle  body,  facilitating  modeling  of  the  anode  heat  transfer.  Flush- 
mounted  planar  single  Langmuir  probes  are  used  for  internal  diagnostics  of  the  boundary  layer 
flow  inside  the  arcjet  nozzle.  Two  nozzle  probes  and  one  constrictor  probe  are  implemented  to 
obtain  ne,  Te,  current  density  and  sheath  voltage  data  in  the  arcjet  for  1 1  <  P/m  <  20  and  flow 
rates  of  60  and  85  mg/sec.  Comparison  between  the  experimental  results  and  a  computational 
model  is  made.  Based  on  the  data  obtained  with  a  nozzle  probe  the  sheath  voltage  is  predominantly 
negative,  so  that  an  electric  field  exists  near  the  anode  which  retards  the  electrons.  The  results 
indicate  a  high  degree  of  ohmic  heating  at  the  sheath  edge,  especially  for  positive  sheath  voltages. 
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The  primary  interest  of  the  Air  Force  for  arcjet  propulsion  is  with  storable  propellants.  The 
research  reported  here  emphasizes  gaseous  mixtures  of  N2  and  H2  for  the  simulation  of  N2H4  and 
H2  arcjets.  Because  the  arcjet  is  a  low  Reynolds  number  device  (Re<1000),  viscous  behavior 
dominates  energy  transfer  processes  in  the  thermal  and  viscous  nozzle  boundary  layers,  which 
couples  with  multi-dimensional  heat  conduction  in  and  from  the  massive  conducting  anode  wall. 
Therefore,  energy  transfer  must  be  described  by  coupled  Navier-Stokes  conservation  relations 
which  depend  in  a  complex  way  on  the  fluid  properties.  Previously,  no  single  computational 
model  could  be  said  to  describe  accurately  the  entire  physics  of  arcjet  flows. 

The  energy  transport  from  the  plasma  located  in  the  equilibrium  constrictor  region  of  the 
arcjet  is  a  strong  function  of  the  number  density  and  distribution  of  electrons  from  ionized  species. 
The  degree  of  non-equilibrium  can  significantly  alter  the  expected  electron  number  density.  The 
diffusion  of  species,  particularly  electrons,  is  likewise  important.  The  electric  current  distribution, 
not  known,  is  strongly  coupled  to  the  ne  and  Tg  distribution  in  the  vicinity  of  the  relatively  cold 
constrictor  and  anode  walls.  Thermal  and  chemical  non-equilibrium  plays  a  dominant  role, 
resulting  in  spatial  variation  in  the  current  attachment  and  anode  heating.  The  proper  description  of 
this  aspect  of  arcjet  nozzle  flow  has  proven  particularly  difficult  to  achieve  in  numerical  models  of 
the  arcjet. 
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I.  A  Plasmadynamics  Model  for  Nonequilibrium  Processes  in  N2/H2  Arcjets 

(AIAA  Paper  No.  95-1961) 


Abstract 

A  general,  two-dimensional,  steady,  non-LTE  (local  thermodynamic  equilibrium)  model  is 
presented  for  nitrogen/hydrogen  arcjet  thrusters.  Both  thermal  and  chemical  nonequilibrium 
processes  are  incorporated  to  describe  a  multi-temperature,  multi-species  N2/H2  plasma.  The 
current  and  anode  temperature  distributions  are  fully  coupled  to  the  plasma  flowfield  with 
realistic  boundary  conditions.  Thermal  nonequilibrium,  chemical  equilibrium  results  are 
presented  for  simulated  hydrazine.  Results  for  both  thermal  and  chemical  nonequilibrium 
hydrazine  and  hydrogen  arcjets  show  significant  thermal  nonequilibrium  near  the  anode,  and 
high  frozen  flow  losses.  It  is  concluded  that  nonequilibrium  chemistry  models  are  necessary  to 
capture  accurately  the  interior  distributions  of  the  plasma  flowfield. 

Nomenclature 

Cj  mean  thermal  speed  (m/s) 

C  coulomb  charge 

Di  effective  diffusion  coefficient  (m^/s) 

elst  elastic  energy  transfer  (W/m^) 

go  gravitational  acceleration  at  sea  level  (m/s^) 

h  Planck's  constant  (J-s) 

hi  enthalpy  of  species  i(J/kg) 

H  heavy  species  total  particle  energy  (J) 

I  integrated  total  current  (A) 

Isp  specific  impulse  (s) 

j  current  density  (A/m^) 
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kb  Boltzmann's  constant  (J/K) 

kc  Coulomb  constant  (Nm^/C^) 

kf  Forward  rate  coefficient  (m^/s  or  m^/s) 
kr  Reverse  rate  coefficient  (m^/s  or  m^/s) 
Le  Electron  inelastic  loss  rate  (W/m^) 

Lh  Heavy  species  inelastic  loss  rate  (W/m^) 
mi  mass  of  species  i  (kg) 

mij  reduced  mass  (kg) 

m  mass  flow  rate  (kg/s) 

M  Mach  number 

ni  number  density  of  ith  species  (m"3) 
p  pressure  (N/m^) 

Pq  integrated  ohmic  power  deposition  (W) 

Qij  collision  cross  section  (m^) 

Ta  anode  temperature  (K) 

Te  electron  temperature  (K) 

Tg  heavy  species  temperature  (K) 

T  thrust  (N) 

u  axial  velocity  (m/s) 

V  radial  velocity  (m/s) 

V  velocity  vector  (m/s) 

V  voltage  (volts) 

w  azimuthal  swirl  velocity  (m/s) 

x,r,6  cyhndrical  coordinates 

Xi  mole  fraction  of  species  i 

yi  mass  fraction  of  species  i 

Xcon  constrictor  length  (mm) 


7 


Xmol  mixture  parameter  in  XmolN2  +  H2 

Zj  thermodynamic  partition  functions 

ai  catalytic  efficiency  of  species  i 

Ed  dissociation  energy  (J) 

Ei  ionization  energy  (J) 

£  anode  emissivity 

Xq.  anode  thermal  conductivity  (W/m-K) 

electron  thermal  conductivity  (W/m-K) 

X.h  heavy  species  thermal  conductivity  (W/m-K) 
Ae  Spitzer  logarithm  term 

|J,  gas  viscosity  (kg/m-s) 

V  collision  frequency  (s‘l) 

6vH2  H2  vibrational  temperature  (K) 

0vN2  N2  vibrational  temperature  (K) 

p  mass  density  (kg/m^) 

CT  electrical  conductivity  (mhos/m) 

X  stress  tensor  (N/m^) 

collision  integral  (m^) 
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I.  Introduction 


Electrothermal  arcjets  offer  a  significant  advantage  in  specific  impulse  and  cost  over 
conventional  satellite  propulsion  systems.  Low  power  1-2  kW  systems  have  been  flight 
qualified  and  are  now  used  for  north-south  station  keeping,  while  higher  power  20-30  kW 
designs  are  in  advanced  development  for  orbit  transfer  applications.  ^ 

The  arcjet  is  shown  schematically  in  Fig.  1.  The  fluid  dynamics  and  energy  transfer  within 
the  arcjet  thruster  are  nonlinear  and  strongly  coupled.  The  propellant  is  injected  upstream  of  the 
constrictor,  typically  with  an  azimuthal  (or  'swirl')  velocity  component.  The  energy  conversion 
mechanism  for  an  arcjet  is  an  arc  discharge,  produced  by  a  voltage  difference  between  the 
thruster  nozzle,  which  is  the  anode,  and  a  conical  cathode  on  the  upstream  side  of  the  constrictor. 
Typical  operating  voltages  for  1  kW  class  arcjets  are  100  V  DC  at  a  current  of  I  =  10  A.  The  arc 
current  distribution  depends  on  several  factors  in  addition  to  the  geometry.  The  distribution  is 
coupled  to  the  propellant  mass  flow  rate,  composition,  thermal  properties,  electrical  conductivity, 
and  gas-dynamic  properties. 

Both  chemical  and  thermal  processes  in  the  arcjet  are  described  as  nonequilibrium.  Arc 
current  is  converted  to  electron  thermal  energy  through  ohmic  dissipation.  The  electrons 
transfer  thermal  energy  to  heavy  species  in  the  arc  plasma  through  collisions.  This  energy  is  then 
converted  to  thrust  as  the  fluid  accelerates  through  the  nozzle.  In  regions  of  low  pressure  and/or 
elevated  ohmic  heating,  collisional  coupling  between  electrons  and  heavy  species  may  not 
establish  equal  gas  and  electron  temperatures.  Additionally,  flow  velocities  are  large  so  that  fluid 
residence  times  in  the  nozzle  are  of  the  order  of  a  few  |xs;  much  of  the  energy  invested  in 
dissociation  and  ionization  is  frozen.  Pressures  drop  from  roughly  1  atm  at  the  constrictor  to 
millitorr  conditions  at  the  exit;  radial  diffusion  of  electrons  from  the  arc  core  is  important  in 
determining  the  arc  stmcture. 

The  thrust  and  specific  impulse  produced  by  an  arcjet  are  determined  by  several  factors 
including  the  power  transferred  to  the  propellant,  the  extent  of  both  kinetic  and  chemical 
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nonequilibrium,  the  propellant  gas  mixture,  and  the  arcjet  nozzle  geometry.  Other  factors  which 
contribute  to  the  overall  performance  include  the  thermal  loading  of  the  anode,^  and  the  voltage 
sheaths  at  the  electrodes.  A  comprehensive  model  of  an  arcjet  thruster  is  highly  complex, 
including  plasmadynamic,  fluid  dynamic,  radiative,  and  surface  and  volumetric  heat  transfer 
phenomena. 

Most  arcjet  design  strategies  to  date  have  been  empirical,  with  design  improvements  based  on 
experimental  observations.  The  need  for  a  more  complete  description  of  the  complex  physical 
processes  has  resulted  in  a  variety  of  numerical  efforts  ranging  from  simplified  one-dimensional 
models^’^  to  more  comprehensive  two-dimensional  descriptions  of  the  gas  dynamics  and  energy 
transfer  processes.9-22  Comprehensive  arcjet  models  are  required  to  (1)  interpret  experimental 
results,  (2)  understand  the  physical  processes  in  regions  of  the  thruster  where  diagnostic 
techniques  are  challenging,  (3)  improve  arcjet  performance  and  (4)  scale  the  devices  to  higher  or 
lower  power  levels. 

11.  Nonequilibrium  Model 

An  axisymmetric,  steady,  laminar,  continuum  flow,  two-temperature  kinetic  nonequilibrium 
model  is  formulated  for  a  direct  current  arcjet  with  flow  swirl,  a  variable  nozzle  geometry,  and  a 
variable  mixture  ratio  of  nitrogen  and  hydrogen.  Flow  swirl  is  modeled  by  including  an 
azimuthal  momentum  equation.  Two  methods  are  used  to  determine  species  populations.  If 
chemical  equilibrium  is  assumed,  the  species  populations  are  computed  from  a  generalized  law  of 
mass  action.  For  chemical  nonequilibrium  simulations,  populations  are  determined  from 
species  continuity  equations. 

The  model  predicts  thrust,  specific  impulse,  and  internal  fields  for  pressure  p,  mass  density 
p,  seven  species  densities  ni,  and  anode  Ta,  electron  Te  and  heavy  species  Tg  temperatures. 
Additionally,  the  model  predicts  the  current  density  distribution  j,  voltage  potential  V,  and 
velocity  components  u,  v,  and  w. 
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The  details  of  the  model  are  presented  in  six  sections  below  as:  (1)  the  fluid  dynamic 
equations,  (2)  the  energy  equations,  (3)  the  electromagnetic  equations,  (4)  the  species  number 
density  equations,  (5)  transport  coefficients,  and  (6)  boundary  conditions.  The  arcjet  geometry 
is  described  as  (x,r,0)  where  x  is  the  axial  coordinate,  r  is  the  radial  coordinate,  and  0  is  the 
azimuthal  coordinate. 


1.  Fluid  Dynamic  Equations 

The  viscous  fluid  dynamic  equations  are  summarized  below: 


Axial  momentum: 


a(puu)  1  a(rpuv)  _4  3/  9u\  i  3/,,  Bui 
3x  r  dr  “sS^axI 

2_5(^  +  ilinr— 1  ^ 

■3S^3r)‘3r  3x  r^^^'axj'ax 


(1) 


Radial  momentum: 

3(puv)  I  3(rpvv)  _  ^  .2v?|£.4HI 

d\  r  3r  3)^  3x  /  3r  ^  3r  |  3r  /  3r  3x  3r  3r  3  ^2 

2  pw2 

3r  ^^^3x  /  3r  r 


Azimuthal  momentum  (3p/30  =  0): 

3(pu(wr))  1  3 d  (  3(wr)\ 

_!:_  +  i^(rpv  +  2n)(wr)) 


(3) 


Axisymmetric  continuity: 

|(pu)  +  ^rpv)  =  0 


(4) 


2.  Energy  Equations 

The  assumption  of  kinetic  nonequilibrium  gives  separate  energy  equations  for  the  electrons 
and  heavy  species.  For  the  chemical  equilibrium  model,  the  electron  energy  equation  is: 
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^  ■((■'e'’  -  i)(CpeTe ))  =  ^  MtJ  +  V  .(^Cp,Tp,Vy.)  -  elst  -  radiation  (5) 

where  'elst'  represents  energy  lost  through  elastic  collisions  with  the  heavy  species, and 
'radiation'  denotes  optically  thin  radiation  loss  due  to  continuum  bremsstrahlung .24  The  elastic 
energy  transfer  source  term  is  calculated  as  the  sum  of  electron-molecule,  -atom,  and  -ion 
contributions  as: 


elst=3k„(T,-Tg)n,m,l;  ^ 

i  e  T 


(6) 


where  the  average  collision  frequencies  v^j  between  electrons  and  heavy  species  are  calculated 
using  the  mean  electron  thermal  speed  and  collision  cross  sections  Qes  as:24 


V  .  =  '\l  n  Q_- 


(7) 


where  the  electron-ion  cross  sections  are: 


Qe. 


ion 


Tte^k^lnAg 


(8) 


Electron-molecule  collision  frequencies  are  multiplied  by  an  energy  loss  factor  5  to  account  for 
inelastic  losses  to  the  internal  energy  modes. 25  a  constant  loss  factor  of  6  =  3000,  chosen  for 
the  simulations  presented  here,  is  discussed  below. 

The  chemical  equilibrium  energy  equation  for  heavy  species  (molecules,  atoms,  and  ions)  is: 

V  KvH)  =  V  -(X^VT  )  +  Z  i^V  (pCCpiTg  +  h?)  D^Vy,)  +  elst  +  V  (v  t)  (9) 

i  e  T 

where 
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(10) 


H  = 


1  e 


H: 


C  T 
pi  g 


+  + 


The  ionization  and  dissociation  energies  of  molecules  and  atoms  are  arbitrarily  included  in  the 
heavy  species  energy  equation  as  reference  enthalpies  hi®.  The  reference  enthalpies  and  specific 
heats  are  summarized  in  Table  1. 

For  chemical  nonequilibrium  simulations,  a  finite  rate  chemistry  model  is  used  to  determine 
inelastic  energy  exchange  rates  between  the  particle  classes.  This  model  is  summarized  in  Table 
2,  which  lists  the  individual  reaction  rates  and  energy  exchanges  assumed  for  each  chemical 
process. 26-29  j^e  species  continuity  equations  are  multiplied  by  the  reference  enthalpies 
associated  with  each  chemical  reaction,  and  the  terms  appearing  in  the  energy  equations  as  the 
convection  and  diffusion  of  ionization  and  dissociation  potentials  are  replaced  with  energy 
transfer  rates.  The  electron  energy  equation  for  chemical  nonequilibrium  is  then 


Iv  -  |)(Cp,T, ))  =  ?  ■(X.VtJ  +  V  .(^C^T,D,Vy,) 


+  —  -  elst  -  radiation  -  T  L. 
a  ® 


(11) 


where  individual  energy  loss  rates  Le  are  summarized  in  Table  2. 
The  chemical  nonequilibrium  heavy  species  energy  equation  is: 


V  .KvH)  =  V  KvtJ+  X  iv  .(pCp,Tj  DiVy,) 


where  the  total  energy  is  now: 


1  ^  e 


+  elst-H  V  -(v  •  ij-XLh 


m;V-v 


(12) 


(13) 
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The  anode  temperature  distribution  is  determined  from  the  energy  equation  for  heat 


conduction  in  a  solid.  Neglecting  ohmic  dissipation  in  the  anode  gives 


(14) 


where  A,a  =  X,a(Ta).^® 


3.  Species  Equations 

For  chemical  equilibrium  (CE)  simulations,  the  species  number  densities  are  given  as  a 
function  of  the  pressure  p,  the  heavy  species  temperature  Tg,  and  the  electron  temperature  Te.  A 
seven-species  plasma  composition  of  N2,  H2,  N,  H,  N+,  H+,  and  electrons  is  assumed. 
Dissociation  and  ionization  equations  are  constructed  from  a  generalized  law  of  mass  action  as 
given  by  entropy  maximization. 31  This  yields  the  following  equations  for  the  nitrogen  and 
hydrogen  species: 


n: 


"H2'  VbTg!  h2  , 


H 


3/2 


,ex 


^H2,ex^H2,rot^H2,vibJ 


n; 


%2 


Hd,N2] 

^,ex  1 

1  ^b^'^g  1 

\  h2  j 
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(15b) 
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1  h2  1 

■,ex 
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In  the  above  expressions,  the  excitation  partition  functions  Zex  are  functions  of  Te,  while 
Zrot  and  Zvib  are  functions  of  Tg.  Details  for  the  calculation  of  the  various  partition  functions 
can  be  found  in  References  24  and  3 1-34. 

The  plasma  equation  of  state,  quasi-neutrality,  and  the  definition  of  the  mixture  parameter 
xmol  provide  the  additional  relations  required  to  solve  for  the  seven  unknown  species.  Note  that 
the  initial  mixture  composition  is  written  as  XmolN2  +  H2  so  that  pure  hydrogen  (xmol  =0)> 
simulated  hydrazine  (xmol  =  1/2)  and  simulated  ammonia  (xmol  =  1/3)  can  be  easily 
investigated: 


_  %  '*'%+ 


-  %+  +  %+ 

(17) 

(18) 

For  chemical  nonequilibrium  (CNE)  simulations,  Eqs.  (15)  and  (16)  are  replaced  by 
individual  species  continuity  equations,  which  are  written  in  terms  of  the  species  mass  fractions 
yi  =  mini/p: 


V  -(pvyi)  =  V  -(pDiVyj)  +  pj  (19) 

Mixture-averaged  mass  diffusion  coefficients  Di  are  employed,  rather  than  multi-component 
diffusion  coefficients.  With  this  method,  Eq.  (19)  is  solved  for  yN2.  yH>  yN.  yN+»  and  ym- 
The  electron  population  is  then  calculated  from  Eq.  (17),  and  yH2  is  determined  from  overall 
mass  continuity,  or  equivalently  Syi  =  1.  The  species  productions  rates,  pj  =  m|hj,  are 
determined  as  a  function  of  Tg  and  Te  from  the  finite-rate  chemistry  processes  summarized  in 
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Table  2.  The  references  for  the  individual  reaction  rates  are  summarized  in  the  right  column  of 
Table  2.  For  two-way  reactions,  the  reverse  reaction  rates  are  extracted  from  the  forward 
reaction  rate  and  the  'equilibrium  constants'  given  by  Eqs.  (15a-d). 

4.  Electromagnetic  Equations 

The  calculation  of  the  electrical  current  distribution  is  critical  to  understanding  arc  physics 
and  poses  a  challenging  problem  in  arcjet  analysis.  The  current  distribution  is  largely 
determined  by  the  flow-field  electrical  conductivity  a,  which  is  a  strong  function  of  the  electron 
number  density  and  temperature  distribution.  Models  which  assume  local  thermodynamic 
equilibrium  (LTE)  must  artificially  elevate  a  in  the  cool  boundary  layer  region  to  avoid  excessive 
voltage  drops  near  the  anode.^T^^  The  artificial  restriction  that  Te  =  Tg  results  in  an 
underprediction  in  the  ionization  (and  therefore  electron  number  density  and  a)  in  the  boundary 
layer  regions  and  a  subsequent  overprediction  in  power  deposition.  Hence,  the  'conductivity 
floor'  is  employed  elsewhere^0,l4  provide  agreement  between  the  experimentally-observed 
and  predicted  power  depositions. 

Diffusion  of  electrons  from  the  arc  core  towards  the  anode  wall  also  elevates  the  electron 
number  density  and  plays  an  important  role  in  determining  the  boundary  layer  electrical 
conductivity.  This  has  been  demonstrated  by  Butler,  Kull,  and  King,^^  who  show  that  the 
incorporation  of  mass  diffusion  provides  a  reasonable  current  distribution  without  imposing 
artificial  restrictions  on  a  or  the  current  attachment  region.  However,  the  two-temperature 
solutions  of  Miller  and  Martinez-Sanchez^^  and  Keefer  et  al.l^  demonstrate  that  significant 
kinetic  nonequihbrium,  with  elevated  Te,  exists  outside  of  the  arc  core  in  the  constrictor  region. 
These  results  indicate  that  both  diffusion  and  increased  ionization  due  to  elevated  Te  are 
important  factors  which  determine  the  boundary  layer  electrical  conductivity. 

The  correct  approach  to  the  arc  attachment  problem  is  to  model  the  anode  as  an  equipotential 
surface  and  allow  the  current  distribution  to  be  a  model  output  which  is  independent  of  artificial 
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restrictions.  This  is  a  self-consistent  approach  to  solving  for  the  current  distribution  and  electric 
field  simultaneously,  for  which  the  current  distribution  is  properly  coupled  to  the  flow-field. 

Assuming  induced  magnetic  fields  are  negligible,  and  neglecting  electron  pressure  gradient 
driving  terms.  Ohm's  law  and  current  conservation  are  used  to  solve  for  j  and  V:25 


I 

<1 

< 

II 

(20) 

Vj  =  0 

(21) 

5.  Transport  Coefficients 

Transport  coefficients  for  this  model  are  calculated  using  mean  free  path  mixture  rules. 

The  viscosity  is  computed  as: 

n: 

X  HjMjj 

J=1 

(22) 

where  Mij  is  defined  as 

rn - 

V  nii  -(2,2) 

(23) 

and  the  pure  species  viscosity  is 

|ij  =  l/2miniCilj 

(24) 

where  Ci  is  the  mean  thermal  speed,  li  is  the  mean  free  path,  and  Q|j 

is  the  energy-averaged 

collision  integral  for  interaction  between  species  i  and  j. 
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The  translational  contribution  to  the  heavy  species  thermal  conductivity  is  calculated  using  an 
equation  similar  to  (24)  with  the  pure  species  conductivity  given  as: 


9Yi-5  \ 

15Yi-15) 


(25) 


The  electron  thermal  conductivity  is: 


n^C  1 


9Ye-5  ] 
15Ye-15| 


N 

I 

j=l 


n-M.. 


(26) 


and  electrical  conductivity  is  calculated  from: 


a  = 


Uee- 


me  I 

i  ^  e 


(27) 


The  effective  diffusion  coefficients  Dj  are  calculated  from: 


'  4/i®ij 


(28) 


where  the  Dij  are  the  binary  mixture  diffusion  coefficients  given  as 


36 


D„  = 


/2k|,T(nij+mj)y 


16nQjj 


1/2 


;(1,1)\  Jcmjmj  j 


(29) 
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The  electron  and  ion  diffusion  coefficients  are  modified  to  account  for  the  ambipolar  electric 
field.  Assuming  ~  ambipolar  diffusion  coefficients  for 


H+ 


N+ 


electrons,  H+,  and  N+  are: 


£)  -  Pe(M-H+’^H+  M-N+%+)  M-e(PH+%+  Pn+%+) 
M^e^e  M'H+%+  +  J^N+^N+ 


£)  _  +  M'N+^h-)  +  M^H+(Pe^e  '  ^N+^N+) 

M'e^e  +  M-H+^h-  +  M'N+%+ 


D, 


a,N+  -  ■ 


PN^M'e^^e  +  P-H+^H+)  +  M-N+(Pe^e  ~  Ph+%+) 
P^e^e  Ph+"h+  ■'■  Pn+%+ 


(30) 


(31) 


(32) 


where  the  mobilities  |J,  are  calculated  from  |j,j  = 

Vi 

The  above  relations  require  the  energy-averaged  collision  integrals  for  28  species 
interactions.  Collision  integrals  and  collision  cross  sections  are  obtained  from  References  37-43. 
Temperature-dependent  curve  fits  are  used  for  the  transport  property  and  elastic  transfer 
calculations. 


6.  Boundary  Conditions 

The  physical  boundary  conditions  for  the  continuity,  momentum,  energy,  and 
electromagnetic  equations  are  summarized  in  Table  3  and  Fig.  1.  At  the  inflow  boundary  lA,  the 
flow  is  axial  and  the  total  temperature  Tt  and  total  pressure  pt  are  specified.  The  total  pressure  is 
adjusted  to  obtain  the  desired  mass  flow  rate.  The  u- velocities  and  the  inlet  static  temperatures 
are  then  computed  from  the  isentropic  stagnation  property  relations.  The  effects  of  flow  swirl 
are  incorporated  by  specifying  the  azimuthal  inlet  w- velocities  as  a  fraction  of  the  u-velocities. 
At  the  solid  cathode  1C  and  anode  ID  surfaces,  no-slip  conditions  are  employed.  Zero  radial 
gradients  are  specified  on  the  center-line  IE. 
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At  the  outflow  boundary  IB,  the  static  pressure  is  extrapolated  from  the  interior  of  the  flow. 
The  outflow  u-velocities  are  updated  to  satisfy  continuity  using  velocity  corrections  as  given  by 
discretized  forms  of  the  momentum  and  continuity  equations  for  the  exit  plane.  This  is 
consistent  with  the  numerical  solution  method  discussed  in  section  El  below.  Note  that  a  zero 
electron  temperature  gradient  is  employed  at  the  walls  1C  and  ID.  At  low  temperatures,  the 
electron  thermal  conductivity  approaches  zero  so  that  this  boundary  condition  will  have  little 
effect  on  the  results.44 

The  boundary  conditions  for  the  voltage  are  shown  in  Table  3.  The  anode  is  modeled  as  an 
equipotential  surface,  while  the  cathode  voltage  is  equipotential  with  the  voltage  updated  to 
obtain  the  specified  total  current  I.  Note  that  in  this  formulation  the  electrode  sheath  voltage 
drops  Vs  have  been  uncoupled  from  the  bulk  plasma.  The  axial  gradient  of  the  voltage  is  set  to 
zero  at  the  inflow  and  outflow  boundaries  so  that  current  is  forced  to  attach  within  the  arcjet 
nozzle.  Recent  radial  electric  field  and  a  measurements  along  the  exit  plane,  showing  the 
existence  of  roughly  1  AJcvcfi  current  density,  suggest  that  this  boundary  condition  may  need 
refmement.'^^ 

The  boundary  conditions  for  the  anode  energy  equation  are  also  shown  in  Table  3.  On  the 
outer  surfaces  2B  and  2D,  the  local  heat  flux  is  given  by  radiation  to  the  surroundings  at  a 
specified  background  temperature  Tbg  «  Tg  and  a  constant  emissivity  of  e  =  0.31.  Along  the 
anode  surface  2C,  the  net  heat  flux  is  calculated  as  the  sum  of  (1)  the  heat  flux  due  to  conduction 
and  species  diffusion  from  the  flow-field,  (2)  the  radiation  loss  through  the  exit  plane  to  the 
surroundings,  and  (3)  the  sheath  losses  given  by  the  product  of  the  local  current  density  and  an 
assumed  sheath  voltage  drop  Vs  =  13  V.  Note  that  radiation  exchange  between  the  plasma  and 
the  anode  is  neglected  in  this  formulation.  The  upstream  anode  surface  2A  is  either  a  fixed 
temperature  or  zero  heat  flux  boundary. 

The  chemical  nonequilibrium  model  requires  boundary  conditions  for  the  species 
populations.  At  the  anode  surface,  the  diffusion  velocities  are  equated  to  a  fraction  ai  of  the 
species  thermal  velocities:^^ 
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(33) 


‘an  4  V  Tcnij 

The  anode  is  assumed  to  be  catalytic  for  recombination,  so  that  am-  =  ocn-i-  =  an  =  (Xn  =  1- 
Assuming  that  molecules  are  reflected  from  the  surface  gives  aH2  =  ttN2  =  0-  The  inflow  is  in 
chemical  equilibrium  up  to  the  constrictor  entrance,  where  the  chemical  nonequilibrium 
simulation  begins.  The  mass  fractions  at  the  exit  plane  are  extrapolated  from  the  interior  of  the 
flow. 

III.  Numerical  Solution  Method 

The  solution  method  is  structured  around  the  global  continuity  and  momentum  equations. 
These  are  solved  using  a  compressible  form  of  the  pressure-based  PISO  algorithm,^^  where  the 
density  variations  are  implicitly  included  in  the  pressure-correction  procedure  as  described  by 
Rhie.*^”^  Staggered  velocity  and  scalar  grids  yield  a  well-connected  pressure  field.  Upstream 
weighted  densities  are  employed  to  insure  numerical  stability .^8  The  equations  are  discretized 
over  finite  cell  volumes,  and  power  law  differencing  of  convective  and  diffusive  portions  of  the 
linearized  coefficients  preserves  numerical  stability  in  regions  of  the  flow  where  convection 
dominates.  The  solution  process  is  iterative,  and  values  for  the  field  variables  from  the  latest 
iteration  are  used  to  linearize  the  source  terms.  Under-relaxation  is  employed  so  that  updated 
solutions  are  fractionally  added  to  the  solution  from  the  previous  iteration.  The  governing 
equations  are  transformed  into  natural  coordinates  and  solved  on  a  uniform  computational  mesh. 

Figure  2  illustrates  the  solution  algorithm.  The  energy  equations  are  solved  to  update 
gas,  electron,  and  anode  temperature  fields.  The  species  number  densities,  transport 
coefficients,  and  source  terms  are  then  updated  using  the  latest  available  temperatures  and 
pressures.  The  continuity  and  momentum  equations  are  solved  to  update  the  pressure  and 
velocity  fields.  Finally,  the  voltage  and  current  density  distributions  are  updated.  Iterations  are 


21 


continued  until  (1)  the  normalized  energy  equation  residuals  are  within  1%  of  the  electrical  power 
input,  (2)  the  difference  between  the  inlet  and  exit  plane  mass  flow  rates  is  less  than  1%  of  the 
total  mass  flow  rate,  and  (3)  the  axial  and  radial  momentum  equation  residuals  are  within  1%  of 
their  respective  exit  momentum  fluxes. 
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IV.  Results  and  Discussion 


Overview  of  Previous  Simulations 

Before  proceeding  to  the  present  results,  a  discussion  of  the  effect  of  the  electron-molecule 
inelastic  loss  factor  5  on  our  initial  voltage  predictions  is  appropriate.  For  the  NASA-Lewis  1 
kW  arcjet  geometry,  we  found  that  the  6-factors  quoted  in  Ref.  25  (6  «  0(10-100))  resulted  in 
current  attachment  in  the  constrictor,  with  corresponding  underpredictions  of  the  operating 
voltage  of  V  »  40  V  at  I  =  10  A.  We  also  noted  that  increasing  the  mass  flow  rate  had  little  or  no 
effect  on  the  attachment  location.  We  electrically  insulated  the  constrictor  and  observed  little 
effect  on  the  voltage  prediction.  Martinez-Sanchez  recently  re-examined  the  5-factors  for  N2, 
and  values  as  high  as  6  «  4000  were  noted.^  We  choose  6  =  3000  because  (1)  reasonable 
convective  effects  result,  with  the  arc  being  swept  further  downstream  as  the  mass  flow  rate  is 
increased,  as  was  originally  postulated  in  Ref.  8,  and  (2)  better  agreement  is  achieved  between 
the  predicted  and  experimentally  observed  operating  voltage,  while  maintaining  the  equipotential 
anode  voltage  boundary  condition  in  the  constrictor. 

The  constrictor  has  been  electrically  insulated  by  others  to  provide  better  agreement  between 
measured  and  predicted  voltage.  These  models  simulated  the  German  TTl  thruster,  a 
higher  power  10  kW  arcjet  with  an  =  5  mm  long  constrictor,  and  predicted  that  nearly  all  of  the 
electriccd  power  was  deposited  in  the  constrictor.  The  predictions  for  our  geometry,  where  Xcon 
=  0.25  mm,  indicate  that  roughly  50%  of  the  power  is  deposited  downstream  of  the  constrictor. 

Results 

We  present  three  simulations  for  the  geometry  of  a  NASA-Lewis  1  kW  arcjet.  These  are 
summarized  in  Table  4,  where  mass  flow  rate,  voltage,  current,  thrust,  and  specific  impulse  are 
shown.  The  simulations  are  designated  as  either  chemical  equilibrium  (CE)  or  chemical 
nonequilibrium  (CNE).  For  each  case,  the  upstream  anode  surface  2 A  (Fig.  1)  temperature  is 
constant  at  Ta  =  1000  K,  and  the  cathode  1C  temperature  varies  linearly  from  Tc  =  1000  K  at  the 
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inlet  to  Tc  =  3000  K  at  the  cathode  tip.  The  constrictor  diameter  is  0.63  mm.  The  cathode  gap 
spacing,  as  measured  axially  from  the  converging  portion  of  the  anode  surface,  is  Xcat  =  0.58 
mm  for  the  hydrazine  cases  and  Xcat  =1-8  mm  for  the  hydrogen  case.  An  anode  voltage  drop  of 
Vs  =  13.0  V  is  assumed,  and  the  anode  emissivity  is  e  =  0.31. ^2  inlet  flow  swirl  velocity  is 

30%  of  the  axial  inflow  velocity.  The  computational  grids,  shown  in  Fig.  3,  are  comprised  of  62 
axial  by  17  radial  nodes  for  the  plasma,  and  62  axial  by  9  radial  nodes  for  the  anode  domain. 
Simulation  run  times  are  several  hours  on  a  Convex  C-240  mainframe  computer. 

1 .  Hydrazine  Simulations 

The  CE  and  CNE  hydrazine  simulations  are  systematically  compared  in  Figs.  4-15. 
Contours  of  current  density  and  anode  temperature  are  indicated  in  Fig.  4  for  CE  hydrazine.  At 
the  cathode  tip,  the  current  density  is  j  «  40,000  A/cm^.  A  peak  anode  attachment  current 
density  of  j  «  18  A/cm^  is  noted  at  a  location  of  x  «  1 1  mm,  or  roughly  5  mm  downstream  of  the 
constrictor.  The  anode  temperature  increases  from  the  upstream  boundary  value  of  Ta  =  1000  K 
to  Ta  =  1300  K  near  the  exit  plane.  The  CNE  hydrazine  results  are  shown  in  Fig.  5.  The  peak 
anode  current  density  of  j  «  15  A/cm^  occurs  at  x  =  8  mm,  upstream  from  that  for  the  CE  case. 
The  anode  temperature  is  somewhat  lower  than  for  CE,  with  a  value  of  Ta  ~  1200  K  at  the  exit 
plane. 

The  marked  difference  in  the  current  attachment  predictions  is  more  clearly  illustrated  in  Fig. 
6,  where  the  anode  current  densities  are  indicated.  For  CNE,  a  more  diffuse,  bimodal 
distribution  is  noted,  with  the  attachment  closer  to  the  constrictor.  Consequently  the  voltage 
prediction,  indicated  in  Table  4,  is  significantly  lower  at  87  V  for  CNE  versus  1 13  V  for  the  CE 
simulation.  The  anode  Te  distributions  are  indicated  in  Fig.  7.  A  high  degree  of  thermal 
nonequilibrium  is  noted,  with  maximum  electron  temperatures  of  Te  =  8,000  K  and  Te  «  12,000 
K  for  the  CE  and  CNE  simulations  respectively. 

Heavy  species  and  electron  temperature  contours  in  the  constrictor  region  are  indicated  in 
Figs.  8-11.  In  Fig.  8,  a  maximum  Tg  ~  22,000  K  occurs  near  the  center  of  the  constrictor  for 
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CE  hydrazine.  The  central  region  of  the  arc  is  near  thermal  equilibrium,  as  indicated  by  a 
maximum  Te  ==  23,500  K  shown  in  Fig.  9.  The  results  for  CNE  hydrazine  are  similar,  with 
maxima  of  Tg  «  20,000  K  and  Te  =  23,000  K  indicated  in  Figures  10  and  11. 

Axial  velocity  contours  are  shown  in  Figs.  12  and  13.  A  maximum  centerline  velocity  of  u  == 
8  km/s  occurs  at  x  =  8  mm,  and  then  decreases  because  of  heating  and  viscous  losses  to  u  «  6 
km/s  at  the  exit  for  CE  hydrazine.  For  CNE  hydrazine,  a  maximum  velocity  of  u  »  7  km/s  is 
predicted  to  occur  at  x  =  9  mm,  decreasing  to  u  =  6.5  km/s  at  the  exit. 

Finite-rate  chemistry  must  be  included  to  capture  frozen  flow  effects.  In  Figs.  14  and  15, 
exit  plane  predictions  of  electron  number  density  ng  are  presented.  The  predicted  centerline 
density  for  CE  hydrazine  is  ng  =  1.6  x  10^^  cm'^,  while  for  CNE,  the  density  is  3  orders  of 
magnitude  higher  at  Ug  =  1.6  x  10^^  cm"^.  Also  shown  are  Tg  and  Tg  profiles.  More  thermal 
nonequilibrium  is  indicated  for  the  CE  case  with  Tg/Tg  =  2.0  at  the  centerline,  while  the  CNE 
case  is  near  equilibrium  at  Tg/Tg  ~  1 . 1 . 

2.  Hydrogen  Simulation 

The  results  for  the  CNE  hydrogen  simulation  are  indicated  in  Figs.  16  -  20.  Anode 
temperature  and  current  density  contours  are  shown  in  Fig.  16.  The  peak  anode  current  density 
of  j  =  16  Aleve?-  occurs  at  the  attachment  location  of  x  =  8  mm.  Due  to  the  higher  specific  power 
deposition,  the  thermal  loading  of  the  anode  is  somewhat  more  severe  than  for  the  hydrazine 
cases,  with  a  predicted  temperature  near  the  exit  plane  of  Ta  ~  1400  K.  Profiles  of  Tg  and  j  at  the 
anode  surface  are  shown  in  Fig.  17.  The  results  are  qualitatively  similar  to  the  CNE  hydrazine 
case,  with  high  thermal  nonequilibrium  of  Tg  »  12,000  K  predicted  in  the  current  attachment 
region. 

Heavy  species  and  electron  temperature  contours  are  shown  in  Figs.  18  and  19.  In  Fig.  18, 
a  peak  temperature  of  Tg  «  13,500  K  is  noted.  The  electron  temperature  is  similar  in  the  central 
region  of  the  flow  with  Tg  =  14,000  K  indicated  in  Fig.  19.  At  a  location  of  x  =  8  mm,  Tg 
decreases  in  the  radial  direction  reaching  a  minimum  of  Tg  =  3700  K  at  r  =  0.5  mm.  The  electron 
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temperature  then  increases,  eventually  reaching  a  value  of  Tg  «  12,000  K  at  the  anode  surface. 
The  electron  temperature  profile  is  qualitatively  similar  to  the  nonequilibrium  predictions  of 
Miller  and  Martinez-Sanchezl^  and  Keefer  et  al.^^ 

Axial  velocity  contours  are  shown  in  Fig.  20.  A  maximum  centerline  value  of  u  =  15  km/s 
occurs  at  X  =  7  mm  decreasing  to  u  ==  12  km/s  at  the  exit  plane. 

V.  Conclusions 

A  detailed,  propellant-flexible,  N2/H2  arcjet  model  is  developed,  and  is  sufficiently 
comprehensive  to  be  used  for  arcjet  design  studies.  A  thermal  nonequilibrium,  chemical 
equilibrium  model  was  first  developed,  and  then  generalized  to  chemical  nonequilibrium  by 
including  finite  rate  chemistry  and  mass  diffusion.  The  incorporation  of  nonequilibrium 
processes  permits  the  current  distribution  to  be  coupled  to  the  plasma  flowfield  with  realistic 
boundary  conditions  and  without  artificial  modifications  to  the  electrical  conductivity. 

Chemical  equilibrium  and  nonequilibrium  results  have  been  compared  for  hydrazine.  High 
near-anode  thermal  nonequilibrium  is  indicated  in  both  cases  with  Tg/Tg  «  8  and  Tg/Tg  =  12  for 
CE  and  CNE  predictions  respectively.  In  the  constrictor,  the  central  region  of  the  arc  is  near 
thermal  equilibrium  with  Tg  =  Tg  =  20,000  K.  Mass  diffusion  and  finite  rate  chemistry 
dramatically  affect  the  arc  structure.  The  incorporation  of  these  chemical  nonequilibrium 
processes  reduces  the  arc  length  and  lowers  the  predicted  operating  voltage.  Additionally,  higher 
frozen  flow  losses  are  noted.  Predicted  electron  densities  at  the  arcjet  exit  plane  are  three  orders 
of  magnitude  higher  with  nonequilibrium  chemistry. 

Chemical  nonequilibrium  results  have  been  presented  for  hydrogen.  The  results  are 
qualitatively  similar  to  those  for  hydrazine.  However,  lower  arc  temperatures  are  noted,  with 
centerline  Tg  =  Tg  ==  14,000  K  in  the  constrictor. 

The  combined  effects  of  viscous  losses  and  supersonic  heating  were  noted  in  all  cases. 
Centerline  axial  velocities  reach  a  maximum  a  few  millimeters  downstream  of  the  constrictor,  and 
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then  decrease  to  exit  plane  values  of  u  «  6  km/s  and  u  =  12  km/s  for  hydrazine  and  hydrogen 
respectively. 

An  anode  temperature  model  was  included.  For  CE  hydrazine,  maximum  anode 
temperatures  of  Ta  ~  1300  K  are  predicted  near  the  exit  plane.  For  CNE  hydrazine  and 
hydrogen,  maxima  of  Ta  ==  1200  and  Ta  «  1400  are  predicted. 

The  model  has  generated  much  useful  data;  however,  unresolved  issues  remain.  The  voltage 
predictions  largely  depend  on  the  arc  attachment  location.  We  have  found  that  parametrically 
increasing  electron-molecule  inelastic  energy  transfer  increases  the  arc  length.  A  more  general 
physical  model  of  electron-molecule  inelastic  transfer  is  required  to  describe  accurately  the  near¬ 
anode  transport  properties.  Additionally,  a  detailed  electrical  sheath  model  is  needed  to  predict 
near-electrode  voltage  drops. 

We  postulate  that  the  electron  energy  distribution  may  be  non-maxwellian  due  to  the  high 
level  of  inelastic  electron-molecule  energy  loss.  We  are  considering  a  more  complete  description 
of  the  coupling  between  the  electron  energy  and  molecular  internal  energy  modes.  This  may  be 
required  to  describe  the  arc  attachment  and  near-anode  electrical  conductivity  more  accurately. 

The  model  must  be  compared  with  detailed  experimental  measurements  of  plasma  properties 
for  validation  of  both  local  and  global  performance  predictions.  Predictions  from  this  model  are 
compared  with  our  experiments  in  Refs.  49  and  50. 
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VIII.  Figures  and  Tables 
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Table  1:  Species  reference  enthalpies  and  specific  heats. 
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Table  2:  Finite  rate  chemistry  model.  — >  one  way  reactions,  =  two  way  reactions  with  the 

reverse  rate  evaluated  using  the  chemical  equilibrium  constants,  M  is  any  heavy  species  third 
body. 
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Table  3:  Boundary  conditions  for  the  fluid  dynamic,  energy,  and  electromagnetic 


equations.  Energy  equation  boundary  conditions  are  shown  for  both  the  fluid  and  anode 


domains.  For  radiative  boundary  conditions,  Tg  is  the  local  anode  surface  temperature. 
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Table  4:  Summary  of  predictions.  CE  denotes  chemical  equilibrium  and  CNE  denotes 


chemical  non-equilibrium  simulations. 
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Figure  2;  Numerical  solution  algorithm.  The  solution  method  is  iterative,  employing 
successive  under-relaxation  of  the  indicated  field  variables. 


constriaor  region 


6  7 


Figure  3:  Computational  grids.  The  grid  is  tuned  axially  for  high  density  in  the  constrictor  and 
radially  for  high  density  near  the  electrode  surfaces. 
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Cuirent  Density,  A/m^ 


Figure  4;  Anode  temperature  and  current  density  for  hydrazine;  chemical  equilibrium  (CE) 
simulation.  The  maximum  anode  current  density  is  j  «  18  A/cm^  at  x  =  1 1  mm. 
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Figure  5:  Anode  temperature  and  current  density  for  hydrazine;  chemical  nonequilibrium 
(CNE)  simulation.  The  maximum  anode  current  density  is  j  «  15  PJcnfi  at  x  «  8  mm. 
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Figure  6:  Current  density  at  the  anode  surface  for  hydrazine.  Results  are  shown  for  chemical 
equilibrium  (CE)  and  chemical  nonequilibrium  (CNE)  simulations. 
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Figure  7:  Electron  temperature  at  the  anode  surface  for  hydrazine.  Results  are  shown  for 
chemical  equilibrium  (CE)  and  chemical  nonequilibrium  (CNE)  simulations. 
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Figure  8:  Heavy  species  temperature  for  hydrazine;  chemical  equilibrium  (CE)  simulation. 
The  maximum  is  Tg  »  22, (XK)  K  at  the  constrictor  centerline.  The  exit  plane  is  located  at  x  = 
18.5  mm. 
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Figure  10:  Heavy  species  temperature  for  hydrazine;  chemical  nonequilibrium  (CNE) 
simulation.  The  maximum  is  Tg  =  20,000  K  at  the  constrictor  centerline.  The  exit  plane  is 
located  at  x  =  18.5  mm. 


Figure  12:  Axial  velocity  contours  for  hydrazine;  chemical  equilibrium  (CE)  simulation.  The 
maximum  is  u  *  8  km/s  at  x  *  8  mm.  The  centerline  velocity  decreases  to  u  *  6  km/s  at  the  exit 
plane. 
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Figure  13:  Axial  velocity  contours  for  hydrazine;  chemical  nonequilibrium  (CNE)  simulation. 
The  maximum  is  u  =  7  km/s  at  x  =  9  mm.  The  centerline  velocity  decreases  to  u  «  6.3  km/s  at 
the  exit  plane. 


r,  mm 

Figure  14:  Radial  distributions  of  electron  density,  electron  temperature,  and  heavy  species 
temperature  at  the  exit  plane  for  hydrazine;  chemical  equilibrium  (CE)  simulation. 
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Figure  15:  Radial  distributions  of  electron  density,  electron  temperature,  and  heavy  species 
temperature  at  the  exit  plane  for  hydrazine;  chemical  nonequilibrium  (CNE)  simulation. 
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Figure  16:  Anode  temperature  and  current  density  for  hydrogen;  chemical  nonequilibrium 
(CNE)  simulation.  The  maximum  anode  current  density  is  j  «  16  A/cm^  at  x  =  8  mm. 
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Figure  17:  Electron  temperature  and  current  density  at  the  anode  surface  for  hydrogen; 
chemical  nonequilibrium  (CNE)  simulation. 
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Figure  18:  Heavy  species  temperature  for  hydrogen;  chemical  nonequilibrium  (CNE) 
simulation.  The  maximum  is  Tg  =  13,500  K  at  the  centerline  near  the  cathode  tip.  The  exit  plane 
is  located  at  x  =  18.5  mm. 
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Figure  19;  Electron  temperature  for  hydrogen;  chemical  nonequilibrium  (CNE)  simulation. 
The  maximum  is  Te  =  14,000  K  at  the  centerline  near  the  cathode  tip.  High  near-anode  thermal 
nonequilibrium  is  indicated  downstream  of  the  constrictor,  with  Te  *  12,000  K  at  the  anode 
surface.  The  exit  plane  is  located  at  x  =  18.5  mm. 


Figure  20;  Axial  velocity  contours  for  hydrogen;  chemical  nonequilibrium  (CNE)  simulation. 
The  maximum  is  u  =  15  km/s  at  x  =  8  mm.  The  centerline  velocity  decreases  to  u  *  12  km/s  at 
the  exit  plane. 
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n.  Measured  Plasma  Properties  at  the  Exit  Plane  of  a  1  kW  Arcjet 

(AIAA  Paper  No.  95-3066) 

Abstract 

With  recent  advances  made  in  low  power  (1-2  kW)  arcjet  numerical  modeling,  validation 
of  numerical  models  through  experimentation  has  become  increasingly  important.  We  describe 
multiple  electrostatic  probe  techniques  for  characterizing  flow  conditions  at  the  exit  plane  of  a  1 
kW-class  hydrazine  thruster.  An  improved  electrostatic  quadruple  probe  technique,  which  is 
also  modified  to  account  for  plasma  gradients  over  the  finite  probe  volume  and  multiple  ion 
species,  provides  measurements  of  ne,  Te  and  ui/cm  at  the  thruster  exit.  Quadruple  probe  results 
are  also  presented  for  measurements  at  several  centerline  locations  immediately  downstream  of 
the  thruster  exit  plane.  A  spatially-resolved  time-of-flight  electrostatic  probe  technique  is 
employed  for  measurements  of  radial  profiles  of  the  plasma  axial  velocity  ui.  The  quadruple 
probe  results  (Te,  ne,  ui/cm),  coupled  with  independent  measurements  of  plasma  velocity,  are 
used  to  estimate  ion  temperature  Ti=Tg  and  the  extent  of  thermal  non-equilibrium  at  the  thruster 
exit.  Exit  plane  experimental  data  for  ui,  Te  and  Tg  are  compared  with  computational  arcjet 
model  predictions,  showing  agreement  for  ui  and  Tg.  The  model  underpredicts  Te  and  a  direct 
comparison  of  ne  results  is  precluded  by  large  axial  gradients  at  the  thruster  exit. 

Nomenclature 

A  Electrode  geometric  surface  area  [m2] 

Cm  Most  probable  thermal  speed  [m/s] 
e  Electron  charge  [coulombs] 

f  Ion  composition  parameter  (=nH+/ne) 

I  Probe  electrode  current  [A] 

larc  Arqet  operating  current  [A] 

ji  Ion  saturation  current  density  [A/m2] 
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k  Boltzmann  constant  [J/K] 

m  Particle  mass  D^^g] 

ih  Propellant  flow  rate  [mg/s] 

n  Particle  density  [cm'3] 

P  Thruster  power  [kW] 

Pk  Designation  for  electrode  "k" 

Q  Collision  cross  section  [m2] 

rp  Electrode  radius  [mm] 

s  Electrode  clearance  [mm] 

T  Temperature  [K,  eV] 

ui  Ion  flow  velocity  [m/s] 

V  Particle  thermal  speed  [m/s] 

V  Probe  electrode  voltage  [Volts] 

Vd  Electrode  relative  bias  voltage  [Volts] 

V  f  Floating  electrode  voltage  [Volts] 

Vp  Plasma  potential  [Volts] 

Varc  Arcjet  operating  voltage  [Volts] 

X  Axial  distance  from  thruster  exit  [mm] 

xh+,n+  Wake  effect  collection  area  parameters 

5  Flow  divergence  angle  [Degrees] 

£d  Dissociation  energy  [eV] 

r  Mathematical  Gamma  function 

6  Probe  angle  [Degrees] 

K  Correction  to  ne  from  N+  ions 

p.  Equal  to  (mH+/misj+)l/2 

XL  Electrode  "end  effect"  parameter 

X,D  Debye  length  [m] 
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^mfp  Collision  mean  free  path  [m] 
a  Standard  deviation  [mm] 

([)  Equal  to  e/kTe,  [V"l] 

Subscripts 

e,i,n  Electron,  ion,  neutral  species 
g  Gas 

1, 2,3,4  Electrode  designation 

I.  Introduction  and  Background 

Although  research  in  the  area  of  electrothermal  thrusters  is  relatively 
mature,  only  recently  has  the  arcjet  become  applied  to  commercial  satellites.! 
With  these  thrusters  beginning  to  find  a  niche  in  industry,  accurate  numerical 
models  become  increasingly  important  as  design  tools  to  help  in  the 
imderstanding  of  fundamental  arcjet  physics.  As  such,  adequate  measurements  of 
fundamental  plasma  properties  must  be  available  for  the  purpose  of  validating 
numerical  efforts.  It  is  the  purpose  of  this  research  to:  1)  provide  fundamental 
data  for  validation  of  H2/N2  computational  arcjet  models;  2)  add  to  the  existing 
base  of  fundamental  data  in  the  very  near-field  plume  region;  and  3)  demonstrate 
the  feasibility  of  electrostatic  probing  techniques  as  useful  diagnostics  in  the  small 
geometries  associated  with  the  1  kW  arcjet.  In  this  paper  we  present  electrostatic 
time-of-flight  and  quadruple  probes  results  for  the  exit  plane  of  a  1  kW  hydrazine 
arcjet.  Where  applicable,  comparisons  are  made  with  experimental  data  and 
computational  model  predictions  performed  by  our  group. 

Electrostatic  probes  have  found  much  use  in  experimentally  determining 
plasma  parameters  in  space  and  laboratory  plasmas.  Several  authors2-6  have 
employed  classical  Langmuir^  single  probes  in  the  characterization  of  low  power 
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arcjet  plumes.  While  single  probes  are  appealing  due  to  their  relatively  simple 
construction,  interpreting  the  probe  V-I  characteristic  can  be  difficult.  A  viable 
alternative  is  the  triple  electrostatic  probe,®  which  allows  simultaneous 
measurement  of  Te  and  ne  without  the  necessity  of  a  voltage  sweep  and  a 
corresponding  probe  V-I  characteristic.  The  utility  of  triple  probes  for  electric 
propulsion  devices  has  been  demonstrated  in  MPD  thruster  plumes^'l^  and  in  a  1 
kW  hydrogen  arcjet  plume.^®  Additionally,  a  general  review^^  of  electrostatic 
probe  techniques  and  their  uses  in  electric  propulsion  devices  has  recently  been 
presented. 

Under  certain  conditions,  the  crossed-probe  technique^^-lS  can  be  used  to 
measure  ui/cm  (=plasma  velocity /ion  most  probable  thermal  speed)  in  flowing 
plasmas.  The  quadruple  electrostatic  probe^^  combines  the  crossed-probe  and 
triple  probe  techniques  for  simultaneous  measurements  Te,  ne/  and  ui/cm-  The 
utility  of  the  quadruple  probe  has  been  demonstrated  in  an  MPD  thruster^^'^O 
and,  more  recently,  in  the  very  near-field  plume  of  a  1  kW  hydrazine  arcjet. 21 

11.  Quadruple  Probe  Technique 

Figures  la.  and  lb.  show  an  electrical  schematic  and  potential  plot  of  the 
quadruple  probe,  consisting  of  three  cylindrical  electrodes  (Pi,  P2,  P3  )  which  are 
aligned  with  the  plasma  flow  vector  and  one  electrode  (P4)  which  is  perpendicular 
to  the  plasma  flow.  P3  and  P4  are  biased  at  constant  voltages  Vd3  and  Vd4  relative 
to  Pi,  and  P2  assumes  the  floating  potential  Vf  of  the  local  plasma  (l2=0).  The 
circuit  comprised  of  electrodes  1,  3  and  4  is  electrically  floating  such  that  ion 
current  collected  at  P3  and  P4  is  balanced  by  the  electron  current  collected  at  Pi, 
thereby  determining  the  potential  difference  Vd2  [=Vd2(Te)]-  The  three  aligned 
electrodes  1-3  are  similar  to  the  triple  probe  and  yield  simultaneous 
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measurements  of  ne  Te-  The  addition  of  the  perpendicular  electrode  P4  allows 
the  measurement  of  the  parameter  uj/cm/  determined  from  the  ratio  I4/I3  for 
crossed  electrostatic  electrodes  with  Vd3  =  Vd4.  For  the  quadruple  probe,  the 
measured  quantities  are  the  electrode  currents  I3  and  I4  and  the  potential 
difference  between  electrodes  1  and  2,  Vd2-  Currents  I3  and  I4  are  determined  by 
measuring  voltage  drops  across  resistors  (R3  =  R4  =  100  i2).  The  electrode  bias 
voltages  Vd3  and  Vd4  are  provided  by  batteries,  with  Vd3=Vd4=12V. 

It  has  been  shown^  that  probe  contamination  can  have  a  marked  effect  on 
Te  measurements  through  drastic  changes  in  the  measured  quantity  Vd2-  In  the 
present  study,  the  probe  was  cleaned  thoroughly  before  each  data  acquisition 
session  using  ion  bombardment  cleaning.13  This  was  accomplished  by  biasing  all 
four  electrodes  24  V  below  facility  ground  to  attract  ions  to  the  probe.  The 
electrodes  were  then  placed  in  the  arcjet  plasma  stream  for  several  ~  1  s  intervals, 
so  that  the  electrodes  would  glow  and  "burn  off"  the  contaminant  layer. 
Quadruple  probe  sweeps  acquired  before  cleaning  always  yielded  electron 
temperatures  that  were  higher  than  the  post-cleaning  values  by  a  factor  of  ~  2. 

Probe  Construction 

A  schematic  of  the  quadruple  probe  used  in  this  study  is  shown  in  Fig.  2. 
The  probe  is  comprised  of  4  individual  tungsten  wire  electrodes  of  0.25  mm 
diameter  with  an  exposed  length  of  2.5  mm  for  electrodes  1-3  and  2.0  mm  for 
electrode  4.  The  location  of  the  perpendicular  electrode  has  been  modified  from 
that  of  previous  probes,21  to  coincide  with  the  axial  location  of  the  probe  center. 
Each  wire  is  mounted  in  round  single-bore  alumina  (AI2O3)  tubing,  which  is  in 
turn  mounted  in,  and  supported  by,  thin-wall  stainless  steel  tubing.  This  assembly 
is  moimted  in  a  larger  stainless  steel  tube  to  facilitate  connecting  the  probe  tip  to 
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the  probe  positioning  mechanism  discussed  below.  Inside  this  larger  tubing  is  a 
four-bore  alumina  tube  that  electrically  insulates  the  four  probe  leads.  A  high 
temperature  ceramic-based  adhesive22  is  used  to  bond  each  of  the  probe  electrodes 
and  components  in  place.  A  small  quantity  of  adhesive  is  placed  at  the  junction 
between  P4  and  its  single  bore  alumina  tube  to  electrically  insulate  part  of  the 
electrode,  thereby  reducing  the  error  associated  with  determining  its  geometric 
surface  area  (Fig.  2.). 

The  center-to-center  separation  between  electrodes  1  and  2  is  1.4  mm.  This 
spacing,  along  with  electrode  lengths  of  2.5  mm,  allows  an  unaligned  quadruple 
probe  to  handle  flow  divergence  angles  6  ~  30°  without  "shadowing"  of  one 
electrode  by  the  other.  In  an  effort  to  minimize  the  effect  of  the  probe  geometry  on 
the  plasma  flow,  a  gradual  transition  between  the  electrodes  and  the  probe  support 
has  been  implemented.  The  radius  of  the  single-bore  alumina  electrode  supports 
is  ~  0.4  mm,  which  minimizes  perturbations  of  the  plasma  flow  at  the  junction 
between  the  electrodes  and  their  support.  The  exposed  length  of  the  single-bore 
alumina  is  ~  5  mm.  Ceramic  adhesive  is  used  to  fill  the  void,  and  ease  the 
transition,  between  the  electrodes  and  the  stainless  steel  probe  support  tubing. 

Probe  Length  Scales 

For  the  conditions  expected  at  the  exit  plane  probe  location  (Te~0.6  eV,  ne  ~ 
4  X  10^2  cm‘3),  the  Debye  length  is  2.9  x  lO*^  cm  and  the  ratio  of  electrode  radius  to 
Debye  length  rp/X,D=45,  so  that  the  ion  sheath  surrounding  the  electrodes  can  be 
classified  as  thin  with  minimal  error.  For  probe  measurements  5  mm 
downstream,  this  ratio  increases  to  rp/>,D  ~  55.  The  clearance  (s)  between 
electrodes  1  and  3  (Fig.  2)  is  s~l  mm  so  that  s/Xd  >  250  and  the  sheaths  on  adjacent 
electrodes  do  not  interact. 
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To  avoid  strong  sensitivity  of  the  ion  current  to  small  misalignments 
between  the  flow  vector  and  the  electrode  axis,  the  so-called  end  effect  parameter's 
[=(L/uiX,D)(kTe/mi)l/2)]  should  be  greater  than  50.  For  our  plasma  conditions, 
this  parameter  is  conservatively  estimated  at  Xl>250.  The  electrodes  should  also 
be  long  enough  that  the  Bohm  sheath  ion  current  collected  along  the  cylinder 
length  (Ib)  is  much  larger  than  the  current  collected  at  the  probe  end  due  to  the 
convection  of  charged  particles  into  the  probe  tip  (It).  A  first  order  analysis  of  the 
ratio  Ib/It  for  the  conditions  at  the  exit  plane  showed  Ib/It  >  30,  so  that  the  error 
introduced  to  the  collected  ion  current  by  neglecting  the  cylinder  tip  is  less  than 
4%. 


The  radius  rp  of  the  tungsten  electrodes  (Fig.  2)  was  chosen  as  a  compromise 
between  the  thin  sheath  assumption  (rp/X.D  »1)  and  the  assumption  of  free 
molecular  flow  over  the  probe  electrodes  (X,infp/rp  »!)•  For  the  seven  species 
plasma  expected  at  the  exit  plane  (H2,  N2,  H,  N,  H+,  N+,  e),  all  relevant  collision 
mean  free  paths  were  calculated.  The  following  assumptions  are  made  in 
calculating  mean  free  paths:  Te=7000  K,  Tg=2500  K,  rig  =  4  x  cm"3,  and  nneutral 
=  3  X  10f3  cm'3  with  10%  dissociation  at  the  thruster  exit.  The  following  (worst 
case)  mean  free  paths  were  calculated:  A.n-n/rp  >20,  A,i-n/rp  >5,  ^i-i/rp  >30,  X-i-e/rp 
=20,  X.e-n/rp  >20,  and  X.e-e/rp  =30.  Clearly,  the  probe  electrodes  are  operating  in  the 
collisionless  regime.  Additionally,  X-mip/X-D  >200  for  all  collisions,  with  an 
overwhelming  majority  of  the  collisions  having  Xmfp/XD  >1000,  indicating  that 
the  collisionless  sheath  assumption  is  appropriate. 

Quadruple  Probe  Theory 
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Quadruple  probe  analysis  assumes:  1)  The  ion  sheath  surrounding  the 
tungsten  wires  is  thin  compared  to  the  wire  radius  (rp  »  X.d);  2)  Each  probe  wire 
and  sheath  are  collisionless  (A,mfp  »  rp  »  ^d)/  3)  The  electron  energy 

distribution  is  Maxwellian;  4)  Velocity  slip  between  plasma  species  is  negligible; 
and  5)  A  wake  region  of  low  ion  and  electron  density  exists  downstream  of  the 
perpendicular  electrode  P4. 

Derivation  of  quadruple  probe  response  is  similar  to  that  of  the  triple 
probe^  and  has  been  presented  else  where.  "21  jn  our  previous  work^l  we 
presented  some  of  the  revisions  necessary  to  allow  the  application  of  triple  and 
quadruple  probe  theory  to  multicomponent  (N-H)  plasmas.  In  the  present  study, 
the  probe  is  applied  to  a  simulated  hydrazine  plume,  with  multiple  species  of  ions 
present,  whereas  most  of  the  previous  applications  of  triple  and  quadruple  probes 
have  been  for  plasmas  with  one  dominant  species  of  positive  ion.  Single 
Langmuir  probes  have  been  used  in  nitrogen  and  hydrogen-based  arcjet  plumes,^" 

3  but  their  use  is  not  contingent  on  knowledge  of  the  species. 

In  the  present  study,  we  present  a  significant  modification  of  the  quadruple 
probe  theory.  Rather  than  utilizing  an  analytic  expression  for  the  current  collected 
at  electrode  P4, 14  is  related  to  the  current  collected  at  P3  through  the  ratio  of  the 
measured  quantities  I3  and  I4.  The  significance  of  this  modification  is  that  it 
completely  removes  errors  in  Te  and  ne  measurements  that  were  due  to 
assumptions  regarding  effective  current  collection  area  of  P4.  For  a  quadruple 
probe  with  electrode  geometric  surface  areas  A3  (=Ai=A2)  and  A4  and  electrode 
biasing  such  that  Vd3=Vd4,  the  following  expres-sion  can  be  derived: 

^  _  1  +  exp((|)Vd3)  -  2exp[(()(Vd3-Vd2)]  (1) 

(l4/l3){exp[(t)(Vd3-Vd2)]-l} 
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Since  Vd3  is  prescribed  and  Vd2  is  experimentally  determined,  Eq.  (1)  determines 
Tg  =e/k{t)  through  iteration. 


The  ion  saturation  current  density  ji  is  related  to  measured  quantities  by 

(I3/A3)  (I  +  I4/I3)  /, 

^  exp((t)Vd2)  - 1 

Eq.  (2)  is  used  to  determine  ng  from  I3, 14  and  Vd2  by  relating  ji  to  the  electron 
density  through  the  Bohm  sheath  analysis.  For  a  plasma  with  more  than  one 
species  of  collected  ion,  the  Bohm  analysis  must  be  modified  to  account  for  the 
collection  of  each  ion  species  via  the  Bohm  sheath.  For  the  region  of  the 
hydrazine  plume  of  interest,  H'*'  and  N+  are  the  dominant  ions  species,  and  the 
Bohm  expression  becomes;^! 


ji  =  e  (kTe)^/^  exp(- 1) 


I  nH+  I  nN+  ^ 


For  ng=nH'^  +  n]sj’*',  Eqs.(2)  and  (3)  yield 


K  +  Il]exp[T] 

U3/I  I3/  ^2/ _ 

e(kTe)^/2[exp((t)Vd2)-l] 


(3) 


(4) 


In  Eq.  (4),  K  represents  the  effect  of  multiple  ion  species  on  the  probe  response  to 
electron  density: 


where  }X  =  (mH+/inN+)^/^-  For  a  plume  contain-ing  only  H+,  nN+  =  0.0  and  K 
reduces  to  unity,  such  that  Eq.  (4)  reduces  to  single  component  (H+)  form. 


Experimental  observations24,25  show  that  low  nitrogen  ion  densities  are 
found  near  the  exit  plane  in  ammonia  and  hydrazine  arcjets.  It  is  suggested  that 
this  is  a  result  of  charge  exchange  collisions  in  the  expanding  nitrogen-hydrogen 


54 


propellant.  It  can  be  shown  (Fig.  3)  that  equilibrium  number  density  ratios  nN /nH 
and  nN+ZnH'*'  are  equal  for  T  <  7500  K  at  p  =  1  atm.  For  this  temperature  range, 
i^N/nH<  0.4,  due  to  the  preferential  dissociation  of  H2  (ed  =  4.48  eV)  over  N2  (ed  = 
9.61  eV).  The  arcjet  exit  plane  N-H  plasma  is  highly  dissociated  and  partially 
ionized,  with  Te  >  Tg,  and  the  estimated  range  of  gas  temperature  is  Tg  <  6000  K. 
Ionization  by  electron  collision  is  negligible,  and  nitrogen  rapidly  undergoes 
charge-exchange  collisions  with  hydrogen,  the  cross-section  for  which  is  on  the 
same  order  of  magnitude  as  N-H  momentum  transfer  collisions.  Thus:  H  +  N+  <- 
— >  H+  +  N,  and  in  the  absence  of  other  ionization  processes,  the  two  charge- 
exchange  reactions  balance  in  steady  state.  Writing  these  volumetric  collision 
rates  in  the  form  nQv,  with  relative  thermal  speeds  vh,N  =  vH/  and  cross 
sections  QH-N+  =  QN-H+  gives: 

nH+  nH  '' 

which  directly  implies  a  low  nitrogen  ion  density  in  the  exit  plane  (Fig.  3).  Thus, 

the  nitrogen  ion  density  is  effectively  determined  by  the  gas  temperature  Tg  and 
not  the  electron  temperature.  For  Tg  <  6000  K,  un+Zuh'*'  <  0.1  and  K  <  1.07, 
indicating  a  ~7%  increase  of  ne  due  to  the  presence  of  nitrogen  ion  probe  current. 

Crossed  Electrostatic  Probes 

The  use  of  crossed  electrostatic  probes  for  measuring  plasma  velocityl5,17 
and  the  implementation  of  the  crossed  probe  technique  into  the  quadruple  probe 
have  been  presented  elsewhere.19,21  por  a  cylindrical  electrode  oriented  at  an 
angle  0  relative  to  the  plasma  flow,  the  collected  ion  current  is  a  function  of  the 
quantity  (uiZcm)sine,  where  Cm  =  (2kTiZmi)l/2,  For  the  case  of  two  equally  biased 
probes,  with  one  probe  aligned  with  the  flow  vector  and  the  other  one  normal  to 
it,  the  single  species  collected  ion  current  ratio  for  the  thin  sheaths  case  is:!^ 
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(7) 


I4-2A4 
l3  i^A3 


OO 


In  order  to  apply  the  crossed  probe  technique  to  the  multicomponent  (H2-N2) 
plasmas  of  interest  here,  significant  revisions  must  be  made  to  Eq.  (7).  In  the 
original  derivation  of  the  crossed  probe  technique,  a  Maxwellian  velocity 
distribution  is  assumed  for  the  single  component  ion  collection.15  Por  the  present 
study,  this  expression  has  been  rederived  using  a  two-component  form  of  the 
Maxwellian  velocity  distribution.  The  derivation  is  similar  to  the  single 
component  case,  and  results  in  a  contribution  to  the  ratio  I4/I3  from  each  of  the 
collected  ion  species.  For  the  two  species  of  ions  assumed  present,  the  ion 
composition  is  defined  by  f  =  nH+/ne  and  (1-f)  =  nN+Zng  (ne=nH+  +  nN+).  The 
resulting  expression  for  the  measured  current  ratio  is 
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where  the  constant  C  is  given  by 


C=- 


2A4 


V7cA3[f  +  (l-f)^] 

In  Eq.  (8),  the  first  term  in  the  brackets  {  }  is  the  H+  contribution  and  the  second 
term  is  the  N+  contribution.  The  most  probable  thermal  speed  for  N+  (cm,N+)  has 
been  written  in  terms  of  Cm,H+  (since  Cm,N+=  H  Cm,H+)-  Note  that  each  of  the 
terms  is  weighted  by  the  relative  density  of  its  corresponding  ion  [f,  (1-f)]. 


(9) 
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Because  a  wake  is  formed  behind  P4,  a  portion  of  the  electrode  area  A4  does 
not  collect  ions.  The  parameters  xh+  and  x]sj+  are  the  fractions  of  the  electrode 
area  A4  available  for  H+  and  N"*"  collection,  respectively.  Because  the  extent  of  the 
wake  effect  is  determined  by  the  relative  magnitudes  of  the  directed  ion  velocity  ui 
and  thermal  speed  of  the  ions,  xh+  and  xn+  are  not  necessarily  equal.  For 
Tg=Ti=2500  K  and  a  directed  velocity  ui  ~  6500  m/ s,  ui/ Cm,H+  ~  1  and  ui/ Cm,N+  ~  4. 
Several  references  discuss  the  wake  effect  for  high  speed  flows. It  has  been 
postulated  and  experimentally  verified  for  an  argon  plasma  at  Ui/ Cm  ~  1,  that  the 
wake  effect  causes  the  effective  collection  area  to  be  half  of  the  geometric  electrode 
area. Since  ui/Cm,H+  is  expected  to  be  ~  1,  xh+  is  assumed  to  be  1/2  so  that  only 
the  front  half  of  P4  collects  ions.  For  the  nitrogen  ions,  the  wake  effect  reduction 
in  collection  area  is  expected  to  be  more  pronounced.  Since  the  directed  velocity  ui 
is  several  times  Cm,N+/  N"*"  is  assumed  to  be  collected  by  the  projected  area^^  of  the 
perpendicular  electrode  only,  such  that  XN+=l/7t- 

Given  f,  Eq.  (8)  becomes  an  expression  for  ui/cm,H+  versus  the  measured 
current  ratio  I4/I3.  Equation  (8)  is  plotted  in  Fig.  4  for  0  <  f  <  1,  xh+=1/2  and 
x]sj+=l/7c.  Pure  H'*’  collection  is  represented  by  f=1.0,  2H+  +  N"'’  collection  by  f=0.67, 
and  pure  N"*"  collection  by  f=0.0. 

If  an  independent  measure  of  the  directed  plasma  velocity  ui  is  known.  Fig. 
4  can  be  used  to  determine  Cin,H+  [=(2kTi/mH+)l/^]  and,  hence,  the  gas 
temperature  Tg=Ti.  We  describe  below  ui  measurements  which  lead  to  an 
estimate  of  Ti. 

Effect  of  Radial  Gradients  on  Probe  Response 
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Quadruple  probe  measurements  at  the  thruster  exit  plane  show  smooth 
symmetric  profiles  for  ne  and  the  electrode  currents  I3  and  I4,  as  shown  below. 
However,  Vd2  is  known  to  have  a  very  asymmetric  profile,^!  which  is  interpreted 
as  a  result  of  the  steep  radial  gradient  in  ne  off-axis,  and  to  a  lesser  extent,  the 
radial  gradient  in  Te-  These  gradients  are  perpendicular  to  the  quadruple  probe 
electrodes  and  the  arcjet  thrust  axis.  In  the  presence  of  plasma  radial  gradients,  the 
individual  electrode  current  equations^O  are  rewritten  in  terms  of  the  plasma 
parameters  at  each  electrode,  and  are  rearranged  to  solve  explicitly  for  Vci2  and  I3 
in  terms  of  radial  profiles  of  Te  and  ne: 


V2=V2  (Te2,K) 

(10) 

Vi=Vi  (Tel,  Te3/  riel,  ne3, 14/I3/  K) 

(11) 

Vd2=V2-Vi 

(12) 

13=13  (Te3/  ne3/  K) 

(13) 

The  subscripts  on  the  parameters  Te  and  ne  in  Eqs.  (10)-(13)  denote  the  local  value 
of  those  parameters  at  each  of  the  electrodes  (1-3)  of  the  quadruple  probe.  The 
ratio  I4/I3  is  a  measured  quantity  and  K  is  a  function  of  the  plasma  composition 
discussed  above.  Since  the  centerline  experimental  data  are  at  zero  gradient 
conditions,  these  data  can  be  used  with  Eqs.  (1)  and  (4)  to  determine  the  centerline 
values  of  Te  and  ne. 

Using  Eqs.  (10)-(13),  Fig.  5  displays  Vd2  vs.  radial  position  for  the  quadruple 
probe  of  Fig.  2,  for  flat  and  Gaussian  ne  and  Te  profiles.  For  Gaussian  ne  and 
constant  Te,  Vd2  increases  monotonically  across  the  thruster.  For  a  Gaussian  Te 
distribution  broader  than  ne,  the  Vd2  profile  shows  more  asymmetry.  The 
constant  Te  case  can  be  solved  analytically  for  Vd2/  giving: 


In  this  equation  the  effect  of  the  density  gradient  can  be  clearly  seen.  As  the  swept 
probe  starts  outside  the  plume  and  enters  it,  nes/nei  <  1.  At  the  axis  nei  =  ne3, 
and  as  the  probe  leaves  the  plume,  nes/nel  >  1-  The  result  [Eq.  (14)]  is  a 
monotonically  increasing  Vci2  across  the  thruster  face. 

Given  arbitrary  profiles  of  Te(r)  and  ne(r),  we  see  that  Vd2(r)  and  l3(r)  can  be 
predicted.  The  inverse  problem,  of  extracting  ne  and  Te  profiles  from  measured 
Vd2  and  I3,  can  also  be  solved.  Since  Vd2  [Eq.  (12)]  is  a  strong  function  of  Te,  and  I3 
[Eq.  (13)]  is  a  strong  function  of  ne,  there  is  sufficient  uncoupling  that  an  iterative 
routine  is  used  to  find  the  unique  Te(r)  and  ne(r)  profiles  which  reproduce  the 
measured  Vd2  and  I3  profiles.  Starting  with  centerline  values  and  assumed 
Gaussian  profiles,  the  routine  marches  the  probe  position  across  the  thruster  face, 
while  generating  a  new  Te  profile.  This  updated  profile  is  used  with  Eq.  (13)  to 
update  ne(r),  and  the  process  is  repeated  until  both  profiles  converge.  Typically 
five  numerical  iterations  are  required  to  reach  convergence. 

III.  Time-of-Flight  Velocimetiy 

Since  in  addition  to  the  quantities  ne  and  Te  a  quadruple  probe 
measurement  yields  the  quantity  ui/(2kTi/mH+)^/^/  knowledge  of  the  heavy 
particle  velocity  ui  can  be  used  to  determine  the  heavy  particle  temperature  Tg  = 
Ti.  We  have  designed  and  implemented  a  time-of-flight  velocimetry  technique 
for  measuring  spatially  resolved  axial  velocities  in  arcjet  plumes.  This  technique 
is  derived  from  a  previous  Current  Modulation  Velocimetry  (CMV)  method. 

In  our  method,  a  short  duration  current  deficit  is  superimposed  on  the  arcjet 
current,  the  effect  of  which  is  monitored  by  an  electrostatic  time-of-flight  (T-O-F) 
probe  as  it  convects  axially  in  the  thruster  plume. 
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Time-of-Flight  Probe 

The  CMV  technique  can  be  made  spatially  resolved  with  electrostatic 
probes,  thereby  eliminating  the  need  to  Abel-invert  the  data  as  required  by  line-of- 
sight  optical  methods.  Figure  6  shows  the  electrostatic  time-of-flight  sensing  probe 
used  in  this  study.  It  consists  of  two  0.75  mm  diameter  tungsten  wires  which  are 
electrically  insulated  along  most  of  their  length  by  two  1.6  mm  o.d.  alumina  tubes. 
The  exposed  length  of  each  tungsten  electrode  is  0.25  mm  and  the  axial  separation 
between  the  two  is  5  mm  (±  2%).  The  elevation  of  the  probe  is  such  that  the  arcjet 
thrust  axis  is  located  between  the  sensing  electrodes.  The  electrodes  are  biased  24 
V  below  facility  ground  (arcjet  anode)  to  monitor  the  local  electron  saturation 
current  density  [~  ne,  (Te)^/^]  at  each  measurement  location.  Variations  in  probe 
electrode  saturation  currents  are  monitored  as  they  convect  over  the  two  T-O-F 
electrodes,  with  the  temporal  separation  of  these  signals,  combined  with  the  probe 
geometry,  indicating  the  flow  velocity. 

Although  this  probe  is  physically  intrusive,  the  probe  electrodes  (0.75  mm 
dia.)  are  in  the  near-free  molecular  flow  regime  and  are  thus  expected  to  have 
little  effect  on  the  flow.  Additionally,  any  plasma  perturbations  due  to  the 
electrodes  (shocks,  collisional  effects,  etc.)  are  expected  to  have  a  similar  result  at 
each  electrode,  thereby  having  minimal  effect  on  the  convective  time  delay. 

Arcjet  Current  Deficit  Circuitry 

The  arcjet  power  processing  unit  (PPU)  produces  current  with  a  high 
frequency  ripple  (~  15-20  kHz,  ±  8%  about  the  mean)  which  causes  "natural" 
fluctuations  in  ne  and  Te.  These  inherent  fluctuations  are  difficult  to  utilize  with 
the  T-O-F  probe  due  to  their  long  period  (~  55  |xs),  compared  with  the  expected 
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convective  time  delay  of  the  T-O-F  probe  (~  1  ps).  Instead,  fluctuations  are 
artificially  introduced  to  the  arcjet  plasma  flow  by  superimposing  a  short  duration 
(~  a  few  jxs)  current  deficit  pulse  on  the  arcjet  operating  current. 

The  circuit  used  to  modify  the  arc  current  is  similar  to  that  used  in  previous 
CMV  studies26  and  is  shown  in  Fig.  7.  The  current  pulse  is  initiated  by  closing  a 
switch  to  a  fast  R-C  circuit  in  parallel  with  the  arcjet,  causing  a  sudden  decrease  in 
arc  current  while  the  capacitor  charges.  This  short  duration  (~  4  |is)  current  deficit 
results  in  a  "tagging"  of  the  plasma  in  the  arc-heating  region.  Switch  closure  is 
accomplished  with  an  IRF-350  field-effect  transistor  (FET),  which  is  closed  by  the 
gate  pulse  shown  in  Fig.  7.  Originally,  the  R-C  circuit  element  was  replaced  by  a 
single  resistor  (20-30  Q),  but  ringing  on  the  edges  of  the  resulting  square  current 
pulse  obscured  the  desired  T-O-F  probe  signals. 

Spatially  resolved  axial  velocities  are  measured  by  triggering  the  arcjet 
current  deficit  pulse  at  several  different  T-O-F  probe  locations.  Since  the  probe  is 
swept  radially  through  the  plume  during  the  experiment,  it  is  necessary  to  trigger 
the  current  deficit  pulse  repeatably  at  each  desired  radial  location.  Probe  location  is 
measured  by  monitoring  the  voltage  of  a  precision  potentiometer  (Fig.  7)  which 
has  been  calibrated  against  the  arcjet  anode  for  probe  position.  The  position 
(voltage)  of  the  desired  measurement  location  is  set  at  the  voltage  comparator 
circuit,  so  that  a  trigger  signal  is  issued  to  the  FIP  214A  pulse  generator  when  the 
probe  reaches  the  desired  radial  location.  In  turn,  the  pulse  generator  supplies  a 
square  pulse  to  "close"  the  FET  and  initiate  the  arcjet  current  deficit.  Because  1) 
the  gate  pulse  must  be  positive  (relative  to  the  FET  voltage),  and  2)  the  arcjet 
cathode  (and  FET)  operates  at  — 100  V  (anode  grounded),  the  FET  Driver  Circuit 
(Fig.  7)  is  a  necessary  buffer  between  the  FET  and  the  HP  pulse  generator.  Since  the 
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probe  is  swept  relatively  slowly  (~  20  cm/s),  the  probe  location  is  effectively 
"frozen"  during  the  ~  3-4  |is  elapsed  time  between  the  initial  trigger  and  the 
monitoring  of  the  tagged  plasma  at  the  T-O-F  electrodes. 


IV.  Experimental  Apparatus 

Electrostatic  probe  experiments  were  carried  out  at  the  exit  plane  of  a  1  kW 
NASA  arcjet  thruster^^  operating  at  a  flow  rate  of  50  mg/ s  of  2H2  +  N2  to  simulate 
fully  decomposed  hydrazine.  Each  component  of  the  propellant  was  individually 
metered  by  Unit  mass  flow  controllers  and  mixed  in  the  propellant  line.  During 
steady  state  operation,  the  thruster  operated  at  10.0  A  and  112  V  (P/m  =  22.4  MJ/kg). 
Arc  current  was  measured  with  a  Hall  effect  current  transducer.  The  thruster 
exhausted  into  aim  diameter  x  1.5  m  long  vacuum  tank^l  at  ~200  mTorr  during 
arcjet  operation.  Anode  temperatures  were  monitored  with  an  optical  pyrometer, 
reaching  ~  1200  K  during  steady  state  operation. 

The  quadruple  probe  and  T-O-F  probe  are  accurately  positioned  and 
repeatably  swept  through  the  thruster  plume  with  the  probe  mount  system  shown 
in  Fig.  8.  The  complete  assembly  is  fastened  to  a  step  motor  (Fig.  8)  mounted  on  a 
linear  translation  carriage.^l  The  motor  and  probe  assembly  carriage  is  capable  of 
±  5  cm  linear  translation  perpendicular  to  the  arcjet  axis,  allowing  the  probes  to  be 
manually  swept  through  the  plume  at  ~  20  cm/s.  Carriage  and  probe  position  are 
determined  to  0.25  mm  by  a  precision  10  kQ  linear  potentiometer  coupled  to  the 
carriage  manual  drive  shaft,  and  calibrated  against  the  outer  diameter  of  the 
anode.  The  probe  support  arm  is  coupled  directly  to  the  motor  shaft  with  the 
probe  tip  on  the  shaft  axis,  allowing  probe  angle  to  be  varied  (in  0.9°  increments) 
without  changing  the  probe  tip  location. 
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Probe  elevation  is  repeatably  located  to  within  0.5  mm  by  aligning  the  probe 
center  with  a  line  etched  on  the  arcjet  anode.  The  separation  between  the  probe 
tips  and  the  thruster  exit  plane  is  reliably  set  to  within  0.1  mm  with  a  spark  gap 
gauge.  The  arcjet  mount  is  known  to  deflect  toward  the  probe  ~0.7  mm  while 
under  vacuum.  Any  error  associated  with  the  arcjet  exit  plane  position  because  of 
this  is  systematic  and  does  not  affect  the  relative  spacing  of  subsequent  axial 
measurement  positions. 

V.  Experimental  Results  and  Analysis 

Experimental  quadruple  and  T-O-F  probe  results  are  presented  for  the  very 
near-field  plume  of  the  arcjet.  To  quantify  the  extent  of  axial  gradients  in  Te  and 
ne,  quadruple  probe  measurements  are  made  at  several  probe  tip  locations  from  1- 
5  mm  from  the  exit  plane.  Radial  profiles  with  probe  tips  located  1  mm  from  the 
exit  are  presented  for  both  the  quadruple  and  T-O-F  probes. 

Quadruple  Probe  Results 

For  the  quadruple  probe,  Te  and  ne  are  determined  by  measurements  of 
Vd2/l3  and  I4.  Since  these  parameters  are  expected  to  be  influenced  by  the 
inherent  current  ripple  associated  with  the  arcjet  PPU,  a  discussion  of  the  effect  of 
the  ripple  on  the  probe  is  warranted.  Figure  9  shows  the  centerline  ion  current 
density  js  (=  I3/A3)  measured  at  electrode  3  and  the  arc  current  larc  versus  time. 
The  current  measured  at  electrode  3  is  in  phase  with  the  arc  current,  except  for  a 
delay  of  ~  3-4  ps  associated  with  the  time  required  for  the  plasma  to  flow  from  the 
arc-heating  region  near  the  cathode  to  the  quadruple  probe.  Accounting  for 
differences  in  average  velocity  and  measurement  location,  this  delay  is  consistent 
with  previous  triple  probe  results  for  the  hydrogen  arcjet.13 
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The  relationship  between  centerline  measurements  of  Vd2  and  js  is  shown 
in  Fig.  10.  While  js  closely  follows  the  phase  of  the  PPU  current  ripple,  Vd2  is  180° 
out  of  phase  with  33.  Since  js  ~  ne(Te)l/2  and  Vd2  increases  monotonically  with  Te 
[Eq.  (1)],  Figs.  9  and  10  indicate  that,  except  for  the  convective  time  delay,  ne 
increases  and  Te  decreases  with  increases  in  larc  during  a  typical  PPU  ripple  cycle. 
Specifically,  Eqs.  (1)  and  (4)  indicate  that  the  ±  10%  current  ripple  causes  an 
associated  ±  18%  ripple  in  ne  (in  phase)  and  a  ±  5%  ripple  in  Te  (180°  out  of  phase) 
about  their  means.  For  the  quadruple  probe  centered  at  x=2.25  mm  from  the  exit 
plane,  the  mean  centerline  values  are  Te  ~0.6  eV  and  ne  ~3.6  x  10^2  cm"^.  The  Te 
results,  although  out  of  phase  with  larc/  are  in  phase  with  the  arc  voltage  Varc 
(which  has  a  ±  3%  ripple)  due  to  the  negative  impedance  characteristic  of  the 
arcjet.  Thus  the  Te  behavior  is  consistent  with  the  variation  of  the  E-field,  which 
heats  the  electrons.  These  results  are  different  from  previous  hydrogen  arcjet 
data,f3  which  indicate  that  ne  decreases  and  Te  increases  with  increasing  larc-  It  is 
unclear  at  this  time  if  this  is  an  indication  of  a  fundamental  difference  in  the 
operating  physics  of  the  hydrazine  and  hydrogen  arcjets. 

Because  of  the  large  variation  in  the  measured  quantities  Vd2/ 13  arid  I4 
with  the  PPU  current,  the  mean  values  of  these  parameters  are  used  to  determine 
the  mean  values  of  Te,  ne  and  ui/cm,H+  through  Eqs.  (1),  (4)  and  (8).  For  radial 
profile  data,  a  moving  average  smoothing  routine  is  applied  to  the  raw  data  so 
that  the  effects  of  the  PPU  ripple  are  eliminated. 

Te  and  ne  Axial  Profiles 

Previous  ne  data  indicate  that  axial  gradients  over  the  length  of  the 
quadruple  probe  may  be  large.4/28,29  Jq  quantify  the  effects  of  gradients  in  Te  and 
ne  over  the  2.5  mm  quadruple  probe  length,  centerline  data  were  acquired  for 
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several  probe  axial  locations.  Five  probe  tip  locations  from  1  to  5  mm  from  the 
thruster  exit,  in  1  mm  increments,  were  utilized.  Results  for  the  measured 
current  density  js  for  five  centerline  axial  locations  are  shown  in  Fig.  11.  The 
geometry  of  the  three  aligned  electrodes  is  also  shown,  with  its  placement 
corresponding  to  the  probe  location  for  the  first  axial  location.  A  50  %  decrease  in 
js  is  indicated  over  the  5  mm  region  investigated,  which  represents  the  combined 
effect  of  changes  in  Te  and  ne,  since  js  ~  ne(Te)^/^. 

Figure  11,  along  with  the  corresponding  Vd2  data  and  Eqs.  (1)  and  (4),  yields 
the  axial  centerline  profiles  of  Te  and  ne  shown  in  Fig.  12.  As  a  check  of  the 
quadruple  probe  and  to  insure  that  the  presence  of  the  perpendicular  electrode 
does  not  influence  the  Te  and  ne  data,  independent  measurements  were  made 
with  a  triple  probe  at  the  first  axial  location  (x=  2.25  mm).  The  quadruple  and 
triple  probe  Te  and  ne  data  were  consistent  within  5  %,  well  within  the 
experimental  error  associated  with  the  experiment.  This  is  not  a  surprising  result, 
since  in  the  limit  I4/I3  ->  0,  Eqs.  (1)  and  (4)  reduce  to  the  triple  probe  equations  for 
Te  and  ne. 

The  electron  temperature  (Fig.  12)  varies  from  6700  to  3600  K  (~0.6  -  0.3  eV) 
between  2.25<  x  <6.25  mm  from  the  thruster  exit.  These  results  are  roughly 
consistent  with  previous  far-field  plume  results  of  Te  ~0.2-0.3  eV  for  30<x<150 
mm  and  Te  ~  0. 1-0.2  eV  for  x=320  mm. 3/4  Further,  the  present  data  seem  to 
corroborate  previous  estimates^S  of  Te  ^  0.5  eV  in  the  arcjet  nozzle  based  on 
arguments  that  Te  is  on  the  order  of  the  atomic  excitation  temperature. 

Results  for  the  electron  density  axial  profile  are  also  shown  in  Fig.  12.  Over 
the  5  mm  range  of  probe  locations  ne  drops  by  ~  30  %,  with  the  magnitude  of 
dne/dx  increasing  slightly  near  the  thruster  exit.  The  ne  results  of  the  present 
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study,  along  with  those  of  other  researchers,  are  summarized  in  Fig.  13.  While  the 
ne  gradient  based  on  the  present  data  is  large  near  the  exit  (dne/dx  ~  -  0.5  x  10^2 
cm'3/mm  at  x=  2-3  mm),  it  is  small  compared  with  the  axial  gradient  indicated  in 
the  nozzle  interior  and  at  the  exit  plane.  It  thus  appears  likely  that  large  axial 
variations  in  the  ne  gradient  (large  d^ne/dx^)  are  present  in  the  thruster  very  near¬ 
field  plume.  The  present  ne  data  appear  consistent  with  previous  results  with  the 
exception  of  previous  data  for  x=0  and  x=10  mm  (Fig.  13).29  Flowever,  these  data 
represent  upper  limits  placed  on  ne  based  on  Stark  broadening  of  line-of-sight 
(non-Abel  inverted)  emission  spectroscopy  signals,  and  thus  are  consistent  with 
our  probe  results. 

Based  on  the  axial  variations  of  Te  and  ne  reported  above,  it  is  apparent  that 
the  quadruple  probe  results  represent  plasma  conditions  averaged  over  the  axial 
dimension  (2.5  mm)  of  the  probe.  Additionally,  the  very  large  ne  gradients 
presumed  to  exist  immediately  downstream  of  the  thruster  exit  require  extra  care 
in  locating  and  reporting  probe  positions  for  both  optical  and  intrusive 
measurement  techniques. 

Exit  Plane  Radial  Profiles 

Traces  for  the  measured  quadruple  probe  voltage  Vd2  and  current  I3  are 
shown  for  a  typical  exit  plane  radial  sweep  in  Fig.  14.  These  data  have  been 
smoothed  so  that  the  traces  shown  represent  the  mean  signals  (with  the  PPU 
ripple  effects  removed).  As  discussed  above,  the  I3  trace  is  symmetric  about  the 
thruster  centerline,  while  the  Vd2  trace  exhibits  the  asymmetry  attributed  to  radial 
gradients  in  ne  and,  to  a  lesser  extent,  Te  over  the  probe  face.^l  The  modified 
quadruple  probe  theory  [Eqs.  (10)-(13)]  is  used  to  solve  for  the  radial  profiles  Te(r) 
and  ne(r)  resulting  from  the  radial  quadruple  probe  data  Vd2  and  I3  (Fig.  14).  This 
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analysis  indicates  that  the  Te  profile  is  very  flat  near  the  centerline  (r<3  mm)  and 
nearly  uniform  at  Te  ~  0.6  eV.  This  uniform  Te,  along  with  Eq.  (13),  is  used  to 
determine  the  radial  ne  profile  at  x=2.25  mm  from  the  thruster  exit,  shown  in  Fig. 
15.  Uncertainty  in  ne  values  far  from  the  centerline  is  larger  due  to  misalignment 
between  electrode  3  and  the  flow  divergence  angle  5. 

T-O-F  Velocity  Probe  Results 

The  effect  of  the  time-of-flight  current  pulse  on  thruster  operation  is 
demonstrated  in  Fig.  16,  which  shows  the  arcjet  current  larc  and  voltage  Varc 
versus  time  during  a  typical  pulse.  The  current  drops  suddenly  as  the  FET  switch 
closes  (Fig.  7)  and  the  0.1  pF  capacitor  begins  to  charge.  As  the  capacitor  nears  the 
end  of  the  charge  cycle  (~  a  few  1  p,s  R-C  time  constants),  the  current  through  the 
FET  leg  of  the  circuit  decreases  so  that  the  arc  current  resumes  the  normal  PPU 
ripple  after  a  small  positive  ring.  The  width  of  the  arcjet  current  deficit  pulse  is  ~  4 
ps.  The  nature  of  the  behavior  of  Varc  is  not  immediately  clear  from  Fig.  16. 
However,  the  derivative  of  the  arc  current  (dlarc/^ii)  h^s  the  same  qualitative 
behavior  as  Varc/  implying  a  possible  inductive  effect  resulting  from  the  arcjet 
geometry  and/or  from  the  current  pulse  circuitry. 

The  T-O-F  probe  response  is  shown  in  Fig.  17  for  a  typical  centerline  axial 
velocity  measurement.  The  signal  measured  at  the  downstream  electrode  has 
been  increased  by  a  factor  of  6  to  account  for  the  lower  electron  current  (lower  Te 
and  ne)  5  mm  downstream  of  the  upstream  electrode.  Both  the  upstream  and 
downstream  electrodes  show  noise  beginning  at  t  ~  2  ps  that  corresponds  with  the 
action  of  closing  the  FET  switch.  Approximately  3  ps  later,  the  effect  of  the  arcjet 
current  pulse  is  monitored  as  a  decrease  in  the  electron  current  measured  by  the 
upstream  electrode.  After  a  time  delay  for  the  tagged  flow  to  traverse  the  5  mm 
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electrode  separation,  a  similar  decrease  is  noted  at  the  downstream  electrode.  For 
both  of  the  signals,  the  width  of  the  decrease  in  the  measured  electron  current  is  ~ 
4  |is,  approximately  the  same  as  that  of  the  current  deficit  pulse. 

Axial  velocities  are  inferred  from  the  convective  time  delay  demonstrated 
in  Fig.  17  and  the  known  separation  (5.0  mm)  between  the  two  T-O-F  electrodes. 
The  time  delays  are  derived  from  the  probe  response  using  the  temporal 
separation  between  the  minima  (due  to  the  decrease  in  electron  current)  measured 
at  each  electrode.  This  is  accomplished  by  taking  the  derivative  of  each 
(smoothed)  signal  as  shown  in  Fig.  18.  The  difference  in  the  locations  of  the  zero 
crossings  is  the  time  required  for  the  tagged  flow  to  travel  the  5.0  mm  electrode 
separation. 

Results  of  T-O-F  probe  centerline  velocity  measurements  are  shown  in  Fig. 
19  for  data  acquired  during  three  individual  periods  after  thruster  ignition  (cold 
start).  The  data  show  that  the  centerline  axial  velocity  initially  averages  ~  6.0 
km/ s  during  the  period  that  the  thruster  is  warming  up.  After  approximately  5-7 
minutes  of  run  time,  the  axial  velocity  averages  ~  6.5  km/s.  These  data  are 
consistent  with  our  previous  results^l  which  indicate  that  exit  plane  ne 
measurements  reach  steady  state  after  ~  4-8  minutes.  During  the  warm-up  period, 
the  measured  velocities  fluctuate  between  5.8  and  6.3  km/s.  Similarly,  the  steady- 
state  centerline  axial  velocities  range  from  6.3  to  6.7  km/s.  These  variations, 
although  approximately  within  the  experimental  error,  may  be  physical.  Previous 
work  utilizing  CMV  for  a  1  kW  hydrogen  arcjet26  has  shown  similar  velocity 
fluctuations.  In  that  work,  it  was  determined  that  the  phase  of  the  PPU  ripple  was 
not  responsible.  Similarly,  the  arcjet  current  deficit  pulse  was  always  triggered  at 
the  same  phase  in  the  PPU  current  ripple  (Fig.  16)  for  the  data  presented  in  this 
study.  It  is  unclear  what  may  be  causing  this  phenomenon. 
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The  T-O-F  probe  was  also  used  for  off-centerline  measurements  to  generate 
an  axial  velocity  profile  at  the  thruster  exit.  For  these  data,  an  attempt  was  made 
to  align  the  probe  with  the  flow  divergence  angle  6.  It  was  found  that  using  a  fixed 
angle  (0=0°)  probe  produced  a  flat  or  even  inverted  profile  with  a  minima  on  the 
centerline.  An  off  axis  T-O-F  probe  that  is  grossly  misaligned  with  the  flow 
divergence  angle  (0<<6)  will  see  a  "faster"  plasma  at  the  downstream  electrode 
than  at  the  upstream  electrode.  Additionally,  the  electrode  separation  along  the 
streamline  is  less  than  the  geometric  separation  (5.0  mm)  by  the  factor  cos(5).  Both 
of  these  factors  artificially  increase  the  measured  axial  velocities,  with  the  former 
phenomenon  having  a  greater  effect  than  the  latter.  Although  the  plasma 
streamlines  are  not  known  a  priori,  the  error  associated  with  T-O-F  probe 
misalignment  for  off  axis  measurements  can  be  minimized  by  rotating  the  probe 
to  align  it  with  the  expected  flow  divergence  at  a  given  location.  The  flow 
divergence  has  been  shown*^  to  be  quite  large  ~  10  mm  from  the  thruster  exit,  and 
is  expected  to  be  large  at  the  exit  as  well.  For  T-O-F  measurements  at  r  =  0,  1,  2,  3 
and  4  mm,  the  probe  angles  used  were  0°,  10.8°,  18.0°,  25.2°  and  45.0°,  respectively. 
It  was  assumed  that  velocities  measured  in  this  manner  were  along  the  probe 
angle  and  were  converted  to  axial  components. 

Figure  20  shows  the  results  of  axial  velocity  measurements  made  at  the 
thruster  exit  plane,  with  the  upstream  electrode  located  1.0  mm  from  the  exit.  The 
profile  is  fairly  symmetric  and  is  compared  with  a  computational  model  as 
discussed  below. 

Analysis  of  Crossed  Probe  Data 
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Since  a  measurement  of  the  centerline  value  of  the  plasma  velocity  ui  has 
been  made,  quadruple  probe  data  can  be  used  to  estimate  the  heavy  particle  kinetic 
temperature  Tg=Ti.  For  the  quadruple  probe  located  at  x=2.25  mm,  the  measured 
ratio  of  the  ion  currents  collected  by  electrodes  3  and  4  is  I4/I3  =  0.55  to  0.60.  This 
current  ratio,  along  with  Fig.  4  for  f=0.9  (K=1.07),  yields  ui/cm,H+  =  0.95  to  1.05, 
which  indicates  that  the  most  probable  thermal  speed  for  H+  (cin,H+)  is 
approximately  equal  to  the  directed  ion  velocity  ui.  For  the  measured  centerline 
velocity  ui=6.5  km/s,  the  heavy  particle  temperature  ITi=(cm,H+)^n^H+/2]  is  Tg=Ti 
~  2500  K  (0.2  eV).  Therefore,  with  the  previous  results  for  Te  ~  6700  K,  the  extent 
of  thermal  non-equilibrium  on  axis  at  the  thruster  exit  is  Tg/Te  ~  0.4. 

Since  the  perpendicular  electrode  P4  in  reality  measures  I4  based  on  a 
distribution  of  ui  and  Ti  along  its  length,  a  discussion  of  the  errors  introduced  by 
the  probe  size  is  in  order.  Because  the  quadruple  probe  theory  for  determining  Te 
and  ne  uses  the  measured  value  of  I4/I3,  no  error  due  to  the  presence  of  the 
perpendicular  electrode  is  involved  in  determining  Te  and  ne.  This  has  been 
verified  (above)  with  independent  triple  probe  measurements.  In  interpreting  the 
crossed  probe  data,  however,  the  geometry  of  the  probe  becomes  a  factor.  The 
current  I4  is  likely  slightly  underrepresented  due  to  changes  in  ne  and  ui  along  the 
electrode  length.  However,  because  ui/cin,H+  is  proportional  to  l/(Ti)l/2,  a  Ti 
profile  peaked  at  the  centerline  has  the  opposite  effect  on  I4.  Assuming  a  10% 
error  in  I4, 14/13=  1.1(0.57)~0.63.  The  ratio  ui/cin,H+  becomes  ~  1.15,  and  Ti  becomes 
2000  K  (0.17  eV). 

Comparison  of  Computational  Model  and  Data 

A  computational  arcjet  modeling  effort^O'^l  has  been  undertaken  by  our 
group  concurrent  with  the  exit  plane  probe  studies  discussed  above.  The  model  is 
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an  axisymmetric,  seven  species  hydrogen /nitrogen  plasma  code  that  utilizes  a 
PISO  algorithm  to  solve  the  computational  domain  up  to  the  exit  plane.  Separate 
energy  equations  are  formulated  and  solved  for  the  electrons  and  the  heavy 
species.  The  anode  temperature  distribution  is  included,  and  the  plasma  electrical 
conductivity  is  coupled  with  the  plasma  properties,  which  allows  a  self-consistent 
solution  for  the  current  distribution.  The  model  is  capable  of  both  chemical 
equilibrium  and  non-equilibrium  simulations.  A  complete  discussion  of  the 
model  and  assumptions  has  been  presented  elsewhere.30,31 

Comparisons  are  made  between  experimental  results  and  model  predictions 
for  a  1  kW-class  arcjet  operating  on  50  mg/s  of  2H2+N2  to  simulate  fully- 
decomposed  hydrazine  propellant.  The  arc  current  is  10.0  A  for  both  model  and 
experiment.  The  specific  powers  for  the  laboratory  thruster  (112  V,  10  A)  and 
numerical  model  (94  V,  lOA)  are  22.4  and  18.2  MJ/kg,  respectively.  The  model 
specific  power  is  contingent  on  the  empirical  value  of  the  anode  fall  potential 
used.  For  the  model  results  discussed  below,  a  chemical  non-equilibrium 
simulation  is  employed. 

The  centerline  exit  plane  ne  predicted  by  the  model  is  ne=3  x  10^^  cm"3, 
which,  given  the  largely  varying  axial  ne  gradients  discussed  above,  seems 
consistent  with  previous  experimental  data  (Fig.  13).  Because  of  these  large 
gradients,  a  direct  comparison  between  the  quadruple  probe  exit  plane  data  (probe 
tip  =  1  mm  from  exit)  and  the  computational  model  is  difficult. 

Centerline  exit  plane  heavy  particle  temperature  predictions  agree  favorably 
with  experiment  (Tg^model  =  2800  K,  Tg^exp.  ~  2000-2500  K),  but  Te  results  differ  by  a 
factor  of  ~  2  (Te,model  =  3000  K,  Te,exp.  ~  6600  K).  Based  on  experimental  results  of 
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the  present  and  previous  studies,  it  thus  appears  that  the  numerical  model  may  be 
underpredicting  electron  temperatures  at  the  thruster  exit  plane.  Note  that  model 
predictions  of  thrust,  ui,  Isp  and  Tg  are  unaffected,  since  the  relative  electron 
concentration  is  low  (ne/nn  <  0.001).  Internal  nozzle  diagnostics32  performed  by 
our  group  are  expected  to  help  resolve  the  discrepancy  between  experiment  and 
numerical  model  predictions. 

Fig.  20  shows  a  comparison  between  the  predicted  axial  velocity  profile  and 
the  results  of  the  T-O-F  probe  surveys.  Agreement  is  excellent  and  within 
experimental  error  in  most  locations.  The  numerically  predicted  profile  is  slightly 
more  narrow  and  peaked  than  the  measured  profile.  This  is  an  expected  result 
because  the  model  flow  is  constrained  in  the  nozzle  (no-slip  at  wall),  which  differs 
from  the  free  expansion  that  occurs  in  the  plume  in  the  vicinity  of  the  T-O-F 
probe. 

VI.  Summary  and  Conclusions 

Experimental  results  for  two  electrostatic  probe  diagnostic  techniques  are 
presented  for  the  very  near-field  plume  of  a  1  kW  hydrazine  arcjet.  An  improved 
quadruple  probe  theory,  showing  excellent  agreement  with  independent  triple 
probe  surveys,  was  used  to  measure  centerline  axial  profiles  of  Te  and  ne  for 
x=2.25-6.25  mm  from  the  exit  plane.  While  the  gradient  in  ne  is  slight  over  the 
region  investigated,  indications  are  that  the  density  varies  by  1-2  orders  of 
magnitude  within  the  first  few  millimeters  downstream  of  the  exit  plane.  For  the 
axial  position  nearest  the  thruster  exit,  Te=6600  K,  ne=3.6  x  10^2  cm"3  and  Tg=Ti  ~ 
2000-2500  K.  With  the  heavy  particle  temperature  measurement,  the  extent  of 
thermal  non-equilibrium  has  been  quantified  at  the  thruster  exit.  A  spatially 
resolved  technique  for  measuring  plasma  axial  velocities  has  also  been  developed 
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and  implemented.  The  radial  velocity  profile  was  measured  at  the  exit  plane,  with 
a  peak  centerline  velocity  of  6500  m/s.  The  above  results  were  used  to  evaluate 
the  predictions  of  a  numerical  H2/N2  arcjet  model.  Agreement  was  very 
favorable  for  both  the  exit  plane  centerline  Tg  and  radial  ui  profile.  Based  on 
previous  and  current  data,  it  appears  that  the  model  may  underpredict  exit  plane 
Te-  Internal  nozzle  diagnostics  are  expected  to  illuminate  the  causes  of  this 
discrepancy.  Comparison  of  ne  predictions  and  experimental  data  could  not  be 
attempted  due  to  the  large  gradients  in  ne  present  at  the  exit  plane.  Clearly,  there 
are  significant  advances  to  be  made  in  understanding  fundamental  arcjet  physics, 
through  both  experimentation  and  modeling  efforts. 
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Fig.  la.  Quadruple  probe  electrical  schematic.  Measured  quantities  are  Vd2/ 13 
and  I4. 
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Fig.  lb.  Quadruple  probe  electrode  potential  plot  for  Vci3=Vd4- 
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Fig.  3.  Ionic  and  atomic  number  density  ratios  for  simulated  N2H4  at  p=l  atm. 


Fig.  4.  Velocity  ratio  ui/cm,H+  vs.  I4/I3  for  0.0<f<1.0.  For  this  study,  f  (= 
nH+/ne)  ~  0.9,  corresponding  to  K=1.07. 


Fig.  5.  Calculated  quadruple  probe  response  to  axisymmetric  Gaussian  profiles 
of  Te  and  ne  with  standard  deviation  a  (mm). 
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Fig.  6.  Side  view  schematic  of  the  time-of-flight  electrostatic  probe. 


Fig.  7.  Schematic  of  circuit  used  to  trigger  arcjet  current  deficit  pulse  for  T-O-F 
velocimetry. 


Ceramic  TC  NASA  1  kW  Arciet 


Fig.  8.  Schematic  of  the  probe  mount.  Probe  axes  can  be  rotated  in  0.9° 
increments. 
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Fig.  9.  Arcjet  current  and  ion  saturation  current  density  showing  the  effects  of 
the  PPU  ripple.  The  delay  between  larc  and  js  is  ~  3-4  ps. 


Fig.  10.  Vd2  and  js  vs.  time.  Figs.  9  and  10  show  that  ne  increases  and  Te 
decreases  with  increases  in  larc- 


Fig.  11.  Axial  variation  of  the  ion  saturation  current  density.  Horizontal  error 
bars  are  2.5  mm  long,  indicating  the  length  of  electrode  3. 
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Distance  from  Thruster  Exit,  mm 


Fig.  12.  Axial  variation  of  Te  and  ne-  The  ne  gradient  is  ~  -0.5  x  10^2  cm‘3/mm 
nearest  the  thruster.  Experimental  error  in  Te  is  ±  15%. 
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Fig.  13.  Comparison  of  centerline  ne  results  from  the  present  work  and 
previous  studies. 
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Fig.  14.  Smoothed  profiles  of  l3(r)  and  Vci2(r)  measured  by  the  quadruple  probe. 
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Fig.  15.  Measured  rie  radial  profile  for  x=2.25  mm  from  the  thruster  exit  plane. 
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Fig.  16.  Variation  of  larc  arid  Varc  during  the  current  deficit  pulse.  The  pulse 
width  is  much  smaller  that  the  period  of  the  PPU  ripple. 


Fig.  17.  Variation  of  the  T-O-F  probe  signals  during  the  current  deficit  pulse. 
Probe  electrode  signal  delay  is  ~  0.77  ps. 
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Fig.  18.  Derivatives  of  the  smoothed  T-O-F  probe  signals  depicting  zero 
crossings  at  signal  minima  locations. 
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Fig.  19.  Measured  centerline  axial  velocity  vs.  time  for  arcjet  warm-up.  The 
numerical  model  prediction  is  also  shown. 
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Fig.  20.  Experimental  measurements  and  model  predictions  of  the  axial  velocity 
profile  at  the  thruster  exit. 
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in.  Application  of  Internal  Langmuir  Diagnostics  to  Arcjets 

(AIAA  Paper  No.  95-2386) 


Abstract 

A  1.4  kW  arcjet  thruster  has  been  designed  for  internal  probing  of  the  nozzle.  The  arcjet 
employs  a  single-piece  nozzle  body,  facilitating  modeling  of  the  anode  heat  transfer.  Flush- 
mounted  planar  single  Langmuir  probes  are  used  for  internal  diagnostics  of  the  boundary  layer 
flow  inside  the  arcjet  nozzle.  Two  nozzle  probes  and  one  constrictor  probe  are  implemented  to 
obtain  ne,  Te,  current  density  and  sheath  voltage  data  in  the  arcjet  for  1 1  <  P/rh  <20  and  flow 
rates  of  60  and  85  mg/sec.  Comparison  between  the  experimental  results  and  a  computational 
model  is  made.  Based  on  the  data  obtained  with  a  nozzle  probe  the  sheath  voltage  is 
predominantly  negative,  so  that  an  electric  field  exists  near  the  anode  which  retards  the  electrons. 
The  results  indicate  a  high  degree  of  ohmic  heating  at  the  sheath  edge,  especially  for  positive 
sheath  voltages. 

Nomenclature 

Ap  Probe  area  [m^] 

c  Most  probable  speed  [m/s] 

cw  Specific  heat  of  tungsten  [J/(kg-°K)] 
e  Electronic  charge  [Coulomb] 

E  Total  charge  particle  energy  [J] 

Eoo  Electric  field  in  plasma  [V/mm] 

Eavg  Average  resistive  electric  field  between  electrodes  [V/mm] 

Ip  Total  probe  current  [|iA] 

Ie,sat  Electron  saturation  current  [|xA] 
larc  Arcjet  applied  current  [A] 

Ie,Ii  Electron/ion  current  [|xA] 

Isp  Specific  impulse  [sec] 

lo  Current  at  zero  volts  [|iA] 

j  Current  density  [A/cm^] 

je,i  Electron/ion  thermal  current  density  [A/cm^] 

k  Boltzmann  constant  [J/°K] 

m  Propellant  mass  flow  rate  [mg/sec] 
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me,i  Electron,  ion  particle  mass  [kg] 
ne  Electron  number  density  [m'^J 
neoo  Electron  density  at  sheath  edge  [m'^] 
rp  Probe  radius  [pm] 

Rs  Shunt  resistor  [Q] 

Te,i  Electron  /ion  temperature  [°K] 

ATw  Probe  temperature  rise  [°K] 

VcX3  Plasma  potential  [V] 

Vf  Floating  potential  [V] 

Vp  Probe  biasing  potential  [V] 

Greek  Symbols 

r(x)  Flux  of  charge  carriers  [m^-sec]"! 

KW  Tungsten  thermal  conductivity  [W/m-K] 

X  Mean  free  path  [pm] 

X.D  Debye  length  [pm] 

Collision  frequency  [sec‘1] 

PW  Tungsten  density  [kg/m^] 
a  Electrical  conductivity  [ohm-m]"l 

(t)s  Sheath  voltage  [V] 

I .  Introduction 

The  goal  of  this  research  is  to  develop  an  understanding  of  arc  attachment  and  electrode 
heating  processes  at  the  anode  of  a  low  power  hydrazine  arcjet.  The  Langmuir  probe  is 
employed  as  a  diagnostic  for  internal  probing  of  the  arcjet  boundary  layer.  Single  planar  flush- 
mounted  probes  are  situated  in  the  constrictor  and  anode  walls  of  a  1 .4  kW  arcjet  thruster. 

Ultimately,  the  high  performance  limit  of  an  arcjet  is  determined  by  anode  heating.  The 
achievement  of  high  levels  of  Isp,  e.g.  above  600  seconds  for  a  hydrazine  arcjet,  requires  P/in 
(power/  propellant  flow  rate)  exceeding  80  MJ/kg  and  results  in  levels  of  anode  heating  which 
can  cause  electrode  failure.  One  of  the  major  contributions  to  anode  heating  is  the  power 
dissipation  in  the  anode  sheath  which  is  proportional  to  the  sheath  voltage.  In  this  paper  a 
diagnostic  probe  technique  to  provide  measurements  of  that  voltage  is  described. 
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Langmuir  probes  have  been  used  extensively  in  shock  tubes  and  in  hypersonic  ionized 
flowsl“4  and  have  been  a  powerful  tool  of  plasma  diagnostics  since  the  early  1920's.^  They  are 
relatively  simple  in  design  and  straightforward  in  implementation  and  can  provide  a  convenient 
method  of  obtaining  local  plasma  properties  ne,  Te,  j  and  sheath  voltage  (t)s,  near  the  anode 
surface.  Such  measurements  will  help  validate  numerical  model  predictions  near  the  anode,  and 
help  determine  whether  these  models  correctly  describe  the  physics  governing  arc  attachment  and 
anode  heating. 

An  understanding  of  arc  attachment  is  an  important  aspect  of  arcjet  models.  However,  due  to 
the  complex  near-anode  physics  and  the  lack  of  experimental  data  (such  as  ne,  Te  and  c)  in  this 
region,  some  models  have  used  heuristic  approaches  toward  arc  attachment,  e.g.  an  artificial 
conductivity  "floor"  near  the  wall  for  current  attachment,6>^  or  an  insulated  constrictor,  forcing 
arc  attachment  in  the  supersonic  region.  ^ 

Investigations  of  plasma  conditions  inside  the  nozzle  and  constrictor  have  been  performed  on 
low  and  medium  power  arcjets.^-l  1  Zube  et  used  drilled  holes  or  a  quartz  window  to 

obtain  spectroscopic  access.  The  radial  profile  of  H2  excitation  temperature,  electron  density  and 
atomic  hydrogen  density  in  the  constrictor  was  measured.  However,  due  to  problems  with  line 
width  determination  near  the  walls  only  the  data  from  the  inner  60-70%  of  the  constrictor 
diameter  gave  reliable  results.  Temperature  restrictions  on  the  window  and  its  effect  on  energy 
transfer  processes  in  the  constrictor  limit  the  utility  of  the  approach. 

Other  internal  diagnostics  work  includes  that  of  Hargus  et  al,^^  who  performed  internal 
emission  spectroscopy  measurements  in  the  nozzle  expansion  region  of  a  26  kW  ammonia  arcjet. 
Three  optical  access  ports  were  equally  spaced  along  the  nozzle  wall.  Atomic  and  ionic  excitation 
temperatures  of  H  and  Nil  were  obtained  as  well  as  electron  number  densities. 

Curran  et  al.  studied  arc  energy  deposition  in  a  segmented  anode  1-2  kW  arcjet  nozzle,  the 
nature  of  arc  attachment,  and  its  effects  on  performance  characteristics  of  the  device.  The  current 
was  found  to  attach  diffusely  to  the  nozzle  wall,  with  more  than  50%  of  the  input  power  to  the 
arcjet  added  in  the  diverging  section  of  the  nozzle.  The  current  distribution  was  found  to  be 
dependent  on  the  mass  flow  rate. 

In  addition  to  measuring  plasma  properties  in  the  arcjet  boundary  layer,  this  work  employs 
Langmuir  probes  in  the  nozzle  to  assist  in  determining  azimuthal  symmetry  of  the  current 
distribution  and  in  understanding  the  physics  of  arc  attachment.  Unlike  work  with  segmented 
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anodes,  the  Langmuir  probe  approach  provides  a  minimum  disturbance  to  heat  flux  and  current 
attachment  processes  at  the  anode  surface. 

II.  Experimental  Approach 

A  1.4  kW  thruster  was  constructed  with  slightly  larger  dimensions  than  the  standard  NASA  1 
kW  thruster.  The  constrictor  diameter  was  increased  to  1.0  mm  (from  0.63  mm)  and  the 
constrictor  length  to  1.5  mm  (from  0.25  mm),  with  an  area  ratio  of  225.  By  eliminating  the 
anode  insert  and  implementing  a  single-piece  thruster  body  (Fig.  1),  the  arcjet  can  accommodate 
an  array  of  Langmuir  probes,  as  shown  in  Fig.  2. 

Use  of  a  monolithic  body  also  facilitates  modeling  of  the  anode  heat  transfer.  The  arcjet  is 
fabricated  from  a  tungsten  alloy  (HD-18:  95%  W,  3.5%  Ni,  1.5%  Cu)  possessing  better 
machining  qualities  than  2%  thoriated  tungsten.  The  cathode  is  3.2  mm  diameter  at  the  rear  of 
the  arcjet,  stepping  up  to  4.8  mm  at  the  thruster  head.  The  cathode  gap  is  set  at  1  mm.  The 
converging  cone  half  angle  is  30°,  while  the  diverging  section  is  20°  half  angle. 

The  initial  experimental  tests  used  three  single  probes  inside  the  arcjet.  These  Langmuir 
probes  are  made  of  tungsten  wire,  surrounded  by  an  alumina  insulator,  and  inserted  into  a 
stainless  steel  tube.  The  tungsten  wire  and  alumina  tubing  are  held  in  place  with  high 
temperature  (4000  °F)  zirconia  adhesive.  Tungsten  wire  was  selected  due  to  its  high  melting 
point  (3395  °C)  and  high  work  function  (4.53  V),  giving  low  electron  emission  from  the  probe. 
A  calculation  of  the  temperature  rise  of  the  probes  due  to  a  probe  current  Ip  at  a  bias  voltage  Vp  is 
given  by; 


ATw  = 


2IpVpt' 


.0.5 


Ap(7tpwCwKw) 

where  Ap  is  the  probe  collection  area,  t  is  time,  and  pw,  Cw,  and  Kw  are  the  probe  material 
properties.  Probe  collection  areas  are  0.10  mm^  for  the  nozzle  probes  and  0.022  mm^  for  the 
constrictor  probe. 


For  a  15  volt  probe  bias  at  a  typical  4  mA  current  for  50  msec,  the  maximum  temperature  rise 
ATw  of  the  constrictor  probe  is  30  °K  and  about  10  °K  for  the  nozzle  probes.  Assuming  the 
probes  operate  at  the  wall  temperature,  these  calculated  temperature  rises  indicate  that  probe 
heating  is  not  a  problem. 
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The  constrictor  contains  one  0.18  mm  diameter  probe  and  the  anode  contains  two  0.38  mm 
diameter  probes  (A  and  B)  separated  by  180°  at  an  axial  location  x  =  2.2  mm  from  the 
constrictor  probe  (Fig.  2). 

Each  probe  is  fastened  to  the  grounded  anode  with  a  miniature  tubing  connector,  and  installed 
in  the  arcjet  body,  facilitating  removal  of  the  probes  for  inspection  or  replacement.  The  tungsten 
probe  wires  are  connected  to  a  subminiature  thermocouple  mount.  Belden  twisted  pair  shielding 
wires  connect  the  mount  to  a  vacuum  connector,  then  to  a  circuit  board,  and  terminate  at  the 
differential  amplifier  of  a  10  MHz  digital  oscilloscope  via  BNC  cables  (Fig.  3).  The  thruster  is 
seated  on  a  mount  inside  a  1.5  m^  vacuum  tank,  with  a  background  pressure  of  ~  100  mTorr  for 
flow  rates  of  50  mg/sec.  The  simulated  hydrazine  propellant  flow  rate  can  be  varied  between  50 
and  90  mg/sec  and  is  controlled  by  two  Unit  mass  flow  controllers.  Plasma  parameters,  ne,  T  e 
and  j  are  derived  from  the  probe  characteristic  data,  obtained  by  biasing  the  probe  with  a  function 
generator  and  measuring  the  probe  current. 

The  function  generator  supplies  each  probe  with  a  ±14.5  volt  peak-to-peak  sinusoidal  10  Hz 
pulse.  This  gated  pulse  is  triggered  externally  by  a  6  VDC  signal.  The  function  generator 
frequency  is  kept  at  10  Hz  to  avoid  unwanted  voltage  drops  across  the  low  inductance,  100 
current  sensing  resistors  Rs.  The  oscilloscope  collects  4096  samples  at  a  rate  of  32  kHz,  twice 
the  frequency  of  the  16  kHz  arcjet  power  processing  unit. 

A  9.4  |J,F  capacitor  is  used  across  each  shunt  resistor  (Fig.  3)  to  filter  noise  in  the  signals. 
The  probe  voltage  is  determined  from  the  function  generator  output  and  is  corrected  for  the 
potential  drop  across  Rs  during  data  reduction. 

Prior  to  each  experiment  the  probes  are  visually  inspected  for  erosion  and  checked  for 
continuity  in  the  circuit.  The  arcjet  is  initially  started  at  60  mg/sec,  taking  about  10-12  minutes  to 
achieve  a  steady  state  temperature  of  ~  790  °C.  Before  the  system  is  triggered  to  initiate  data 
acquisition,  the  probes  are  cleaned  by  applying  a  0.5  Hz,  14.5  volt  square  wave  pulse,  achieving 
probe  cleaning  by  ion  bombardment.  Once  the  arcjet  has  achieved  steady  state  operation  the 
function  generator  is  then  externally  triggered  and  all  three  probes  are  biased  simultaneously. 
The  data  presented  in  this  paper  is  for  flow  rates  of  60  and  85  mg/sec,  and  1 1  <  P/m  <  20. 

Simple  Probe  Theory 
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The  Langmuir  probe  is  a  metallic  electrode  inserted  into  a  plasma,  biased  at  positive  and 
negative  voltages  Vp  with  respect  to  the  plasma  potential,  Voo.  The  current  Ip  collected  by  the 
probe  is  measured  as  a  function  of  Vp,  generating  a  probe  characteristic  (Fig.  4),  providing 
information  on  the  plasma  properties. 

Region  A  (Fig.  4)  is  called  the  electron  saturation  region,  where  Vp  >Voo  so  that  electrons  are 
attracted  to  the  probe  and  ions  are  repelled.  Region  A  provides  the  magnitude  of  the  electron 
saturation  current,  Ie,sat- 

When  the  probe  voltage  equals  the  plasma  potential,  Voo,  there  is  no  net  voltage  gradient 
around  the  probe  and  the  charged  particles  migrate  to  the  probe  because  of  their  random  thermal 
velocities.  In  the  transition  region  B,  low  energy  electrons  are  repelled  and  ions  are  accelerated 
to  the  probe.  If  the  slope  d(lnIe)/dVp  is  constant  then  the  electron  distribution  is  Maxwellian  and 
the  electron  temperature  is  inversely  proportional  to  the  slope. 


The  equation  governing  the  probe  current  collected  in  the  transition  region  (B,  Fig.  4),  is 
derived  by  calculating  the  flux  of  charged  particles  to  the  probe,  assuming  a  Maxwellian  velocity 
distribution  for  the  electrons  and  a  collisionless  thin  sheath.  The  justification  for  the  thin, 
collisionless  probe  sheath  assumption  is  discussed  later  in  the  paper.  The  derivation  for  the 
probe  current  in  the  transition  region  (Appendix  A),  results  in: 


1/2 


exp 


kT, 


-(Vp-V4 


Equation  (2)  is  independent  of  probe  shape,  i.e.  it  holds  for  planar,  cylindrical,  or  spherical 
probes.  The  electron  temperature  is  found  by  taking  the  natural  logarithm  of  Eq.  (2)  and 
differentiating  with  respect  to  the  probe  voltage  to  get: 


d(ln  le)  _  _e_ 
dVp  kTe 


Once  Te  is  obtained  from  the  transition  region,  the  electron  number  density  is  found  from  the 
electron  saturation  current,  Ie,sat>  obtained  from  the  "knee"  in  the  curve  between  A  and  B,  Fig  4. 
This  point  gives  the  plasma  potential  Voo,  and  Ie,sat-  The  electron  number  density  is  found  by 
setting  Vp=Voo  in  Eq.  (2): 

Ip=Ie,sa,  =en,Apy^  (4) 
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giving: 


(all  units  MKS) 

When  the  probe  potential  equals  the  floating  potential,  Vf  (Fig.  4)  the  flux  of  electrons  equals 
the  flux  of  ions  to  the  probe  so  Ip  =  0.  In  region  C  of  Fig.  2,  almost  all  electrons  are  repelled, 
resulting  in  an  ion  sheath.  The  ion  saturation  current  is  much  less  than  the  electron  saturation 
current  because  me  «  mi  and  Te»Ti. 


Ill.  Experimental  Results 
Probe  Characteristic  Data 

Figure  5  shows  a  typical  Ip  vs  Vp  characteristic  for  the  constrictor  probe.  Figure  6  shows  a 
typical  characteristic  for  nozzle  probe  B.  Prior  to  data  reduction  the  probe  voltage  and  current 
data  are  smoothed  using  5  point  moving  average  smoothing. 

Current  Density 

Current  density  j  is  determined  from  the  characteristic  by  finding  the  current  at  Vp  =  0  and 
dividing  by  the  probe  area.  For  the  conditions  in  Figs.  5  and  6,  P/ih=13.8  MJ/kg,  and  ih=60 
mg/sec,  j  =  0.52  A/cm2  for  nozzle  probe  B  and  j  =  0.06  A/cm^  for  the  constrictor  probe. 

Data  Analysis 

The  electron  temperature  is  found  from  the  slope  of  the  transition  region  of  the  In  le  versus 
Vp  plot  (Fig.  7).  To  obtain  the  electron  current  le  from  the  total  probe  current  Ip,  the  ion  current 
li  at  all  voltages  less  than  Voo  is  predicted  and  subtracted  from  Ip,  using  an  analytic  curve  fit  to 
the  ion  saturation  current. 
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The  theoretical  form  of  the  ion  saturation  current  depends  on  the  length  scale  regime,  i.e. 
whether  the  sheath  is  thin  or  thick  with  respect  to  the  probe  radius,  and  collisional  or  collisionless 
with  respect  to  the  charged  particle  mean  free  paths.  Calculations  presented  below  in  Table  I 
show  that  both  nozzle  probes  have  borderline  collisionless  thin  sheaths,  while  the  constrictor 
probe  has  a  thin  collisional  sheath.  The  data  reduction  process  applies  collisionless  theory  to  all 
probes. 

For  the  thin  sheath  case  the  Laframboise  method  is  used.  This  approach  utilizes  a  plot  of 
nondimensional  ion  current  versus  nondimensional  probe  voltage  (Appendix  B).  With  an 
empirical  fit  to  the  Laframboise  plot  of  nondimensional  ion  current  vs  probe  voltage  for  various 
rp/X,D  values,  an  expression  for  the  theoretical  ion  current  is  found  independent  of  Te,  for  Vp  < 
Woo-V 


The  plasma  potential  Voo  is  found  from  a  graphical  analysis  of  the  probe  characteristic.  Ideally, 
the  plasma  potential  is  the  probe  voltage  at  the  "knee"  in  the  characteristic  on  a  semilog  plot  of  le 
vs  Vp,  indicating  the  end  of  the  transition  region  and  the  start  of  the  electron  saturation  regime. 

To  use  Eq.  (6),  an  estimate  of  ne  is  required.  This  is  done  by  solving  Eq.  (6)  for  ne  and 
choosing  a  point  (Vp,  Ip=Ii)  in  the  ion  saturation  region;  e.g.  (Fig.  6)  Vp  =  -14  V,  li  =  -193  |i.A 
and  Voo  =  0.6  V.  This  data  is  substituted  back  into  Eq.  (6)  for  an  estimate  of  ne-  Using  this  ne 
and  Voo  in  Eq.  (6),  a  theoretical  fit  to  li  is  found  for  all  Vp  <  Voo. 

The  resulting  ion  current  li(Vp)  is  subtracted  from  the  total  measured  current  to  get  le.  Points 
le  <  0  are  discarded  in  order  to  obtain  In(Ie)  vs  Vp.  From  this  plot,  the  slope  of  the  electron 
transition  region  is  found  by  using  a  weighted  least  squares  fit  to  the  data.  If  this  region  is  linear 
(Fig.  7),  the  electrons  have  a  Maxwellian  velocity  distribution,  and  Te  can  be  found  from  Eq. 
(3).  If  the  transition  region  is  not  linear  then  the  Te  obtained  from  the  slope  is  really  an  average 
energy  Eavg  =  3/2kTe.  Once  Te  is  known,  the  electron  number  density  is  found  from  Eq.  (5). 

Results  Varying  P/m 
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Data  for  two  flow  rates  are  presented,  60  and  85  mg/sec,  for  8  A  <  larc  ^  10.4  A.  For  the 
conditions  used  in  these  tests  ne  and  Te  data  were  unavailable  from  the  constrictor  probe  because 
a  probe  characteristic  with  distinct  saturation  regions  could  not  be  obtained  (Fig.  5). 

Figures  8  and  9  are  plots  of  electron  temperature  at  the  anode  wall  as  a  function  of  P/m . 
The  data  is  from  nozzle  probes  A  and  B,  located  about  2.2  mm  downstream  of  the  constrictor 
probe,  in  the  supersonic  portion  of  the  arcjet  nozzle.  There  appear  to  be  no  obvious  trends  with 
P/m  ,  though  the  Te  data  for  probe  A  are  consistently  lower  than  that  for  probe  B.  On  average, 
Te  ~  16,600  °K  ±15%  for  a  flow  rate  of  60  mg/sec,  while  for  85  mg/sec  the  average  Te  is  about 
15,000  °K  ±15%. 

Figures  10  and  1 1  present  the  electron  number  density  measurements  for  both  nozzle  probes. 
For  almost  all  values  of  P/ m  the  ne  values  from  probe  B  are  larger  than  that  from  probe  A.  For 
probe  B  and  a  flow  rate  of  60  mg/sec  ne  varies  between  1  x  lOl"^  m"3  and  2  x  10 1^  m"^  for  14  < 
P/m  <  19.  Similarly,  for  an  85  mg/sec  flow  rate  ne  varies  between  4x101^  m"3  and  1.5  x  lOl^ 
m"3  for  probe  B,  for  1 1  <  P/m  <  15. 

Both  the  electron  temperature  and  the  electron  number  density  data  from  nozzle  probe  A  are 
lower  than  those  from  probe  B.  It  was  discovered,  upon  closer  inspection,  that  probe  B  was 
flush  with  the  wall,  while  probe  A  was  recessed  several  probe  diameters,  reducing  its  data  values 
compared  to  probe  B,  as  shown  in  Figs.  8-11.  This  prevented  simultaneous  measurement  of  j  at 
two  azimuthal  locations,  a  technique  which  will  be  used  in  the  future  to  evaluate  azimuthal 
symmetry  of  the  current  attachment. 

Figures  12  and  13,  for  60  and  85  mg/sec  flow  rates,  show  that  for  nozzle  probe  B  the  largest 
j  value  for  60  mg/sec  is  ~  2.8  A/cm^  for  P/rh  =  17.9  MJ/kg.  The  current  density  has  a 
maximum  of  about  3.7  A/cm^,  for  a  P/in  of  12.5  MJ/kg  and  a  flow  rate  of  85  mg/sec. 

There  is  no  obvious  trend  in  the  j  data  for  85  mg/sec.  For  60  mg/sec,  j  appears  to  increase 
with  P/ in .  For  all  P/m  and  both  propellant  flow  rates,  the  j  measured  by  the  constrictor  probe 
is  generally  significantly  less  than  the  j  measured  by  nozzle  probe  B. 

IV.  Discussion 

The  plasma-wall  interaction  in  the  near-anode  region  of  an  arcjet  presents  a  complex  problem 
in  understanding  the  physics  of  arc  attachment.  This  is  due  to:  (1)  the  non-neutral  sheath  layer; 
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(2)  the  effects  of  the  wall  condition,  i.e.  catalytic  surface  or  emissive;  and  (3)  non-equilibrium 
coupling  of  fluid  dynamics  and  electrodynamics.  Understanding  the  near-anode  physics  is 
important  in  determining  the  transition  between  diffuse  and  constricted  modes  of  current 
transfer.  In  the  diffuse  mode  the  arc  is  dispersed  over  a  large  area,  as  opposed  to  the 
constricted  mode  where  the  arc  is  highly  concentrated.  Operation  of  an  anode  in  a  constricted 
mode  leads  to  large  ohmic  heating  and  erosion,  thus  limiting  the  electrode  lifetime. 

Near  the  anode  surface  there  are  several  coupled  physical  processes:  convection,  diffusion, 
conduction,  and  chemical  kinetics.  The  region  close  to  the  anode  can  be  divided  into  several 
sublayers  where  large  temperature  and  velocity  gradients,  non-neutrality  and  chemical  diffusion 
can  play  important  roles  in  determining  local  plasma  properties. 

To  determine  whether  a  thin  collisionless  sheath  surrounds  the  probes  the  following  length 
scales  are  evaluated,  using  experimental  data:  probe  radius  rp,  Debye  length  and  mean  free 
paths  (mfp)  for  e-H2  and  collisions:  Xe-H2»  Electron-H2  collisions  are 

important  when  electrons  are  collected  by  the  probe,  while  ion-H2  collisions  are  relevant  for  ion 
collection. 

Inside  the  sheath  region,  (<  A.D  from  wall)  charge  neutrality  is  violated  and  electrostatic 
forces  control  the  motion  of  electrons  and  ions.  The  sheath  forms  near  an  electrode  so  as  to 
shield  the  main  body  of  the  plasma  from  the  large  electric  fields  that  exist  within  a  Debye  length 
of  the  wall,2®  and  maintain  charge  neutrality  in  the  plasma. 

Table  I  shows  some  typical  values  for  the  relevant  length  scales  based  on  experimental  data. 

Table  I.  Relevant  length  scales  based  on  experimental  data. 


ne  (m-3) 

Constrictor 

1 

Nozzle 

2  X  1017 

Te(K) 

-4,000 

15,000 

rp  (|i,m) 

85 

190 

Xd  (Itm) 

14 

19 

^e-H2  (|im) 

1 

13 

0.10 

2 
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From  Table  I,  for  the  nozzle  probes,  rp/^D  =  10.  so  the  sheath  is  thin.  The  electron  mfp  is 
on  the  order  of  the  Debye  length,  and  we  use  a  collisionless  approach  for  analyzing  the  electrons. 
The  ion  mfp  is  ~0.1  of  the  Debye  length,  so  that  the  ion  current  is  overstated  by  using  the 
collisionless  assumption. 

The  constrictor  ne  and  Te  (Table  I)  are  estimated  from  the  constrictor  probe  characteristic 
(Fig.  5)  showing  ~  80  |xA  of  electron  saturation  current.  The  calculations  show  that  the 
constrictor  probe  is  collisional  for  both  electrons  and  ions,  but  that  the  sheath  is  thin  (rp/X,D  =  6). 

Sheaths 


A  diagram  of  the  potential  distribution  between  the  cathode  and  anode  of  an  arcjet  is  shown  in 
Fig.  14.  The  anode  sheath  voltage  is  shown  negative  with  respect  to  the  plasma  potential  Vcx5,  so 
that  the  sheath  electric  field  is  towards  the  anode,  causing  electrons  moving  toward  the  anode  to 
be  retarded.  A  negative  sheath  voltage  can  exist  if  there  are  large  pressure  and  electron 
temperature  gradients  near  the  electrode,  so  that  the  current  density,  j,  is  less  than  the  electron 
thermal  current  density,  je- 

The  voltage  in  the  sheath  is  found  by  comparing  the  total  current  density  arriving  at  the  anode 
with  the  electron  thermal  current  density  near  the  anode.  Assuming  the  probe  sheath  is 
collisionless,  the  sheath  voltage  is  derived  first  by  writing  the  absolute  magnitude  of  current 
density  j  as: 


j=je+ji  O) 

where  je  and  ji  are  the  electron  and  ion  thermal  current  densities,  respectively  given  as: 

je  =  (l/4)enece  ;  ji  =  ( l/4)enici  (8) 

where  Ce  and  ci  refer  to  the  electron  and  ion  most  probable  thermal  speeds.  At  the  sheath  edge 
the  electron  density  is  given  by  the  Boltzmann  relation,  so  that  rie  =  =  rig^oo  exp[e(t)s/kTe]. 

Substituting  this  relation  for  ne  and  ni  in  Eq.  (8)  and  using  the  results  in  Eq.  (7)  for  je  and  ji,  an 
equation  for  the  sheath  voltage  is  found: 

engopCe  +  engooCil 

4j  J 
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Typically,  the  ion  thermal  current  density  is  a  few  percent  of  the  electron  thermal  current  density 
so  that  this  term  is  dropped,  resulting  in  Eq.  (lO):^! 


erieooCe 


4j 


(10) 


The  potential  drop  of  the  sheath  can  therefore  be  positive  or  negative,  depending  on  the  relative 
magnitudes  of  j  and  je,  and  can  be  determined  experimentally  given  ne,  Te  and  j  near  the  anode, 
as  extracted  from  a  probe  characteristic. 


The  sheath  voltage  (l)s  from  Eq.  (10)  is  presented  as  a  function  of  P/ ifi  in  Figs.  15  and  16. 
It  is  observed  that  for  almost  all  cases  the  sheath  voltage  is  negative  at  the  nozzle  probe  locations. 


The  experimental  results  for  ne,  Te,  j  and  are  now  be  compared  with  computational  model 
results  for  the  conditions  60  mg/sec  flow  rate,  and  larc  =  8  and  10  A.  The  computational  model 
is  an  axisymmetric,  seven  species  H2/N2  plasma  code,  capable  of  both  thermal  and  chemical 
nonequilibrium  simulations. 22  Comparing  the  experimental  data  with  the  numerical  results  for 
Te  in  Fig.  8,  it  is  seen  that  Te  from  the  model  is  slightly  lower  than  Te  from  the  experiment. 
Closer  agreement  between  the  model  and  our  experiment  occurs  at  the  lower  P/ rh  value  of  14.0 
MJ/kg  for  both  Te  and  ne  data  (Table  E). 


Table  E.  Comparison  of  Nozzle  Probe  B  data  and  Numerical  Model,  rh  =  60  mg/sec. 


fare  ~  8A 


Experiment 

Model 

P/m[MJ/kg] 

13.8 

14.0 

ne  [m-3] 

4.4  X  1017 

8.4x  10l7 

Te[K] 

13,600 

10,300 

j  [A/cm2] 

0.52 

11 

(t)s  [V] 

-1.0 

1.46 
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larc  lOA 


Experiment 

Model 

P/m[MJ/kg] 

17.2 

17.1 

ne  [m-3] 

1.3x  10l7 

l.Ox  10l8 

Te[K] 

20,700 

10,000 

j  [A/cm2] 

0.4 

13.8 

<t)s  [V] 

-0.5 

1.44 

A  major  difference  between  the  experimental  data  and  the  model  is  in  the  current  density,  as 
the  model  shows  current  densities  (Fig.  12)  several  times  the  measured  values.  The  nozzle 
internal  surface  area  is  about  5.5  cm^.  For  larc  =  lOA,  this  means  that  javg  ~  1-8  A/cm^,  which 
has  been  observed  with  probe  B,  Figs.  12  and  13.  Therefore,  our  Langmuir  probe  results  are 
consistent  with  an  even  distribution  of  j  in  the  nozzle,  but  based  on  the  model,  j  is  more  narrowly 
attached.  Fig.  17.  This  difference  can  be  resolved  by  installing  probes  in  different  positions. 

Figure  17  shows  the  model  current  density  distribution  for  a  flow  rate  of  60  mg/sec  and  two 
different  arc  currents.  The  model  predicts  that]  is  not  uniform  and  reaches  a  maximum  value  of 
~14  A/cm^  downstream  of  the  constrictor,  but  about  1.5  mm  upstream  of  the  Langmuir  probe 
location.  Therefore,  it  seems  possible  that  the  probes  are  not  at  the  location  of  maximum  current 
density. 

From  our  experimental  data  for  ne,  Te,  and  j,  a  value  for  the  sheath  voltage  has  been 
calculated.  For  most  of  the  experimental  conditions  (()s  is  negative,  i.e.  the  anode  sheath  electric 
field  retards  electron  motion  to  the  anode.  Fig.  14.  This  implies  less  anode  heating  compared 
with  the  case  >  0. 

For  the  experimental  condition  of  85  mg/sec  flow  rate  and  P/rh  =  14.6  MJ/kg,  the  sheath 
voltage  is  slightly  positive,  =  0.1  V.  The  experimental  probe  characteristic  for  this  positive 
sheath  voltage  is  shown  in  Fig.  18  for  probe  B. 

Comparing  Figs.  6  and  18,  the  plasma  potential  Voo  =  0.6  V  for  (|)s  <  0,  and  Voo  =  -1.8  V  for 
(t)s  >  0.  The  entire  Fig.  18  characteristic  is  shifted  to  the  left  for  the  (l)s  >  0  case. 

The  resistive  electric  field  Eoo  =  j/a  can  now  be  estimated  from  Figs.  6  and  18.  The  electrical 
conductivity  is:22 
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(11) 


a  = 


nge^ 


V 

mg  2^  niCgiZg-i 
iv^e 


where  is  the  energy-averaged  collision  integral.  Based  on  the  mean  free  path  calculations 
in  Section  IV,  the  electron-H2  collision  frequency  is  dominate  over  ion-H2  interactions,  so  that: 


V 

;  =  2^  rii^e^e-i  ^e“H2  ^e-N2 

i^e 


(12) 


Using  constant  collision  cross  sections  for  Qg-H2  ^^*1  Qe-N2->  combined  with  the 
experimental  data  of  Figs.  6  and  18,  a  value  of  a  can  now  be  estimated.  Since  j  is  measured,  an 
estimate  of  Eoo=j/a  is  also  made.  The  results  are  summarized  in  Table  m. 


Table  HI.  Data  for  (l)s  <  0  and  (t)s  >  0. 


ne  [m-3] 

13.8  MJ/kg 
4.4x1017 

14.6  MJ/kg 

6.8x1017 

Te  [°K] 

13,600 

13,900 

j  [A/cm^] 

0.52 

2.2 

a  [Q-m]‘l 

0.15 

0.23 

Eoo[V/mm] 

35 

95 

(|)s  [V] 

-1.0 

O.l 

For  <t)s  <  0,  Eoo  =  35  V/mm  and  for  (j)s  >  0,  Eoo  =  95  V/mm.  These  values  are  larger  than  the 
average  electric  field  in  the  arc  of  Eavg  ~  100  V/4mm  arc  length  =  25  V/mm,  indicating  a  high 
degree  of  ohmic  heating  at  the  sheath  edge,  especially  for  the  case  of  positive  sheath  voltage. 


V.  Conclusion 

We  have  presented  preliminary  experimental  data  from  the  anode  boundary  layer  of  an  arcjet 
and  have  estimated  the  sheath  voltage,  the  electrical  conductivity  and  electric  field  values  at  the 
sheath  edge,  using  single  Langmuir  probes. 
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Though  ne  and  Te  data  were  not  available  from  the  constrictor  probe,  it  was  possible  to 
obtain  j  data,  which  showed  that  under  most  conditions  there  was  significantly  more  current 
downstream  in  the  nozzle  boundary  layer  than  in  the  constrictor. 

Based  on  the  data  obtained  from  a  nozzle  probe  it  appears  that  the  sheath  voltage  is  negative, 
so  that  an  electric  field  exists  near  the  anode  which  retards  the  electrons.  There  may  be  a  critical 
applied  current  density  that  leads  to  a  transition  of  negative  to  positive  sheath  voltage.  What 
governs  the  change  in  sign  of  the  sheath  voltage  is  relevant  to  understanding  the  physics  of  arc 
attachment,  since  the  voltage  drop  in  the  anode  boundary  layer  and  sheath  determines  the 
transition  between  diffuse  and  constricted  modes  of  current  transfer.  Also,  the  sign  of  the  sheath 
voltage  can  be  indicative  of  a  high  degree  of  ohmic  heating  at  the  sheath  edge,  especially  for 
(t)s>0. 


Comparison  between  the  experimental  results  and  a  computational  model  show  approximate 
agreement  for  ne  and  Te.  However,  a  discrepancy  is  found  in  the  current  density  data.  It  is 
possible  that  the  nozzle  probes  are  not  located  in  the  region  of  maximum  arc  attachment;  more 
probes  will  have  to  be  placed  in  the  anode  wall  to  address  the  issue  of  arc  attachment. 
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VIII.  Appendix  A 


Derivation  of  the  Current  in  the  Transition  Region  of  the  Probe  Characteristic 

The  analysis  used  in  this  paper  to  obtain  electron  temperature  and  electron  number  density 
from  the  probe  characteristic  is  based  on  the  assumption  of  a  thin  collisionless  sheath.  The 
derivation  of  the  electron  current  in  the  transition  region  follows  that  of  Swift  and  Schar^^  and  is 
presented  here  for  convenience. 

Assuming  a  collisionless  thin  sheath  surrounds  the  probe,  the  flux  density  of  charged 
particles  reaching  the  probe  is:^^ 

r  =  n(rp)u  (rp)  =  J  f(c)  u  dy  (Al) 

where  dy  represents  the  elementary  volume  in  velocity  space,  f(c)  is  the  velocity  distribution 
function,  n(rp)  is  the  total  charged  particle  concentration  at  the  probe's  surface  and  u(rp) 
represents  the  mean  carrier  velocity  normal  to  the  probe's  surface.  In  the  transition  region  the 
majority  of  the  carriers  are  electrons.  Therefore,  the  current  to  the  probe  is  given  by: 

Ip  =  eTpAp  (A2) 

To  determine  the  probe  current  as  a  function  of  probe  potential,  Tp  needs  to  be  evaluated  in 
terms  of  Vp.  In  Eq.  (Al)  f(c)  is  assumed  to  be  a  Maxwellian  velocity  distribution  for  the 
electrons.  For  a  collisionless  thin  sheath,  with  no  ionizing  reactions,  the  total  energy  E,  of  the 
electrons  is: 23 


E  =  ^  meu2  +  eV{x)  (A3) 

where  only  the  velocity  component  perpendicular  to  the  probe  surface  is  considered  and  V(x)  is 
the  potential  relative  to  the  unperturbed  plasma  at  a  point  x  from  the  probe's  surface,  i.e. 
V(x)  =  Vp  -  V«,.  In  the  unperturbed  plasma  V(x)  =  0  and  E  =  (l/2)meu2. 

Substituting  the  Maxwellian  distribution  for  f(c)  into  Eq.  (Al).: 
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where  ui  is  the  required  velocity  for  the  electron  to  overcome  the  potential  barrier  V(x).  Since 
we  want  the  probe  current  as  a  function  of  probe  potential,  Eq.  (A4)  is  recast  in  terms  of  E,  the 
energy  in  the  unperturbed  plasma  (E  =  (l/2)meu2): 


’’■'H-sr'l 

f- 


exp 


:me_u2 


2kTp 


udu  X 


2kTe  . 


i 


dv  I  exp 


m^w2 


2kT, 


dw 


(A5) 


The  first  integration  is  recast  in  terms  of  E;  therefore  Eq.  (A3)  is  rearranged  to  solve  for  u: 


1/2 


(A6) 


Taking  the  differential  of  u  and  substituting  Eq.  (A6)  into  Eq.  (A4)  an  equation  for  T  is  found  in 
terms  of  E  and  V(x).  Performing  the  integration: 


1/2 

exp 


eV(x) 
kTe  . 


(A7) 


where  r(x)  is  the  flux  of  electrons  crossing  the  sheath  edge.  The  flux  density  of  electrons 
reaching  the  probe  is  found  by  setting  V(x)  =  Vp  -  Voo.  in  Eq.  (A7): 


(A8) 


Therefore,  the  electron  current  in  the  transition  region  of  the  probe  characteristic  is  given  by 
substituting  Eq.  (A8)  into  Eq.  (A2): 


kT, 


-(Vp-V4 


(A9) 
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Collisional  Effects 


Collisions  generally  reduce  the  current  to  the  probe  because  of  the  necessity  for  particles  to 
diffuse  up  to  the  probe  instead  of  arriving  there  by  free  flight .24  The  collisionless  equations 
presented  above  place  an  upper  bound  on  the  collected  current.  To  calculate  the  reduction  in 
the  collected  current,  scattering  collisions  must  be  considered. 

Collisional  sheath  effects  are  important  if  the  electron-neutral  or  ion-neutral  mean  free  paths 
are  much  less  than  the  Debye  length.  If  the  probe  sheath  is  collisional  then  the  above  probe 
theory  requires  modification.  For  example,  the  total  energy  of  an  electron  in  the  sheath  would 
not  be  given  by  Eq.  (A3),  but  would  have  additional  energy  loss  terms  and/or  gain  terms  due  to 
ionization  reactions  or  particle  scattering.  The  species  conservation  and  energy  equations  would 
have  to  be  solved  in  addition  to  Poisson's  equation  in  the  sheath. 


IX.  Appendix  B 
Laframboise  Method 

The  ion  current,  which  must  be  known  to  derive  le,  is  an  empirical  fit  to  a  theoretical 
formulation  of  Laframboise  16  for  collisionless  sheaths,  and  is  presented  here  as  a  convenience  to 
the  reader. 


The  thin-sheath,  transition  to  a  thick  sheath  and  the  thick  sheath  regimes  can  all  be  treated  by 
the  Laframboise  method.  Lafram-boise  solved  Poisson's  equation  numerically,  assuming  a 
Maxwellian  distribution  of  ions  at  Ti  and  electrons  at  Te,  for  various  values  of  rp/^D. 


*  *  ^ 

The  Laframboise  approach  utilizes  a  plot  of  versus  Xp  for  several  values  of  rp/AD-  The 
parameter,  ji ,  is  a  dimensionless  quantity  that  multiplies  the  probe  collection  area  and  accounts 
for  the  Bohm  acceleration  in  the  presheath  region.  Xp  is  a  nondimensional  parameter  that 
indicates  how  far  the  probe  is  biased  into  ion  saturation  compared  to  Te,  given  by;  12 


The  ion  current  is  given  by:  12 


>  _  e{V^  -  Vp) 
P  ■  kTe 


(Bl) 
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To  use  Eq.  (Bl)  both  Ue  and  Te  are  required.  A  simplification  is  made  by  an  empirical  fit  to  the 
Laframboise  plot: 

j‘  =  1.127(Xpf^^ 

With  this  empirical  fit,  an  expression  for  the  theoretical  ion  current  is  found  independent  of  Te: 

Equation  (B3)  can  now  be  used  as  an  analytical  fit  to  the  ion  current  data  and  subtracted  from 
the  total  current  to  obtain  le- 


Fig.  1.  The  1.4  kW  thruster  incorporates  a  single-piece  thruster  body  to  accommodate  internal 
probing  of  the  arcjet  boundary  layer. 


Fig.  2.  Anode  design  with  planar  Langmuir  probes  for  internal  ne,  Te  and  current  density 
measurements.  Nozzle  length  is  19.2  mm;  constrictor  length  is  1.5  mm.  The  nozzle  probes  are 
2.2  mm  downstream  of  the  constrictor  probe,  at  an  area  ratio  of  ~  7.  The  constrictor  probe  is  ~ 
0.5  mm  upstream  of  the  constrictor  exit. 
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Probe  Probe  Constrictor 

B  A  Probe 


Fig.  3.  Experimental  schematic  for  three  single  Langmuir  probes. 


Fig.  4.  Schematic  of  a  typical  probe  characteristic,  showing  the  various  regions  of  interest  that 
provide  information  on  the  plasma  properties. 


Constrictor  Probe  Voltage  V^,  [volts] 

Fig.  5.  Experimental  probe  characteristic  for  the  constrictor  probe  for  conditions  P/iri  =13.8 
MJ/kg,  iri  =  60  mg/sec.  Vp  is  measured  with  respect  to  the  anode  potential,  which  is  zero. 
Probe  collection  area  is  0.022  mm^. 
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Nozzle  Probe  Voltage  V^,  [volts] 

Fig.  6.  Experimental  probe  characteristic  for  the  anode  wall  probe  B  for  conditions  P/rri  =13.8 
MJ/kg,  rri  =  60  mg/sec.  Probe  voltage  is  measured  with  respect  to  the  anode  potential,  which  is 
zero.  The  current  density  j  is  lo/Ap.  Probe  collection  area  is  0. 10  mm^. 


Nozzle  Probe  Voltage  V  ,  [volts] 

P 

Fig.  7.  A  semilog  plot  of  le  versus  Vp  for  P/rri  =13.8  MJ/kg,  rri  =  60  mg/sec.  The  transition 
region,  -3.5  <  V  <  0,  is  linear  over  2  decades,  signifying  a  Maxwellian  electron  distribution.  Voo 
and  Ie,sat  are  obtained  from  the  "knee"  in  the  characteristic. 
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Fig.  8.  Electron  temperature  Te  versus  P/rri  for  60  mg/sec.  Probes  A  and  B  are  located  2.2  mm 
downstream  of  the  constrictor  exit,  and  are  separated  by  180°.  The  values  from  a  numerical 
model22  are  also  shown.  Data  from  probe  A  is  not  valid,  as  discussed. 
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Fig.  9.  Te  versus  P/m  for  85  mg/sec.  Data  from  probe  A  is  not  valid,  as  discussed. 
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Fig.  10.  Electron  number  density  ne  versus  P/ ih  for  a  fixed  flow  rate  of  60  mg/sec.  The  values 
from  a  numerical  model22  are  also  shown.  Data  from  probe  A  is  not  valid,  as  discussed. 
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Fig.  11.  Electron  number  density  ne  versus  P/ih  for  85  mg/sec.  Data  from  probe  A  is  not 
valid,  as  discussed. 
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Fig.  12.  Current  density  at  anode  potential  versus  P/rn  for  60  mg/sec.  The  values  from  a 
numerical  model22  are  also  shown. 
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Fig.  13.  Current  density  at  anode  potential  versus  P/ rh  for  85  mg/sec. 


Fig.  14.  Potential  distribution  diagram  between  the  cathode  and  anode  of  an  arcjet,  and  in  the 
near-anode  region  showing  a  negative  sheath  voltage. 
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Fig.  15.  Sheath  voltage  at  anode  potential  versus  P/ ih  for  85  mg/sec. 
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Fig.  16.  Sheath  voltage  at  anode  potential  versus  P/ih  for  60  mg/sec.  Values  from  a  numerical 
model22  aj-g  also  shown. 


Fig.  17.  Current  density  distribution  along  the  anode  wall  based  on  the  numerical  model22  for 
the  60  mg/sec  and  larc  =  8  and  10  A.  Constrictor  entrance  is  located  at  5.7  mm  and  the 
constrictor  exit  is  at  7.2  mm. 
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Fig.  18.  Probe  B  characteristic  for  (t)s  =  0.11  V  at  85  mg/sec,  P/iri  =  14.6  MJ/kg.  The  location 
of  the  "knee"  makes  Voo  <  0.  Compare  this  characteristic  with  Fig.  6,  for  (|)s  =  -1.0  V. 
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